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Abstract 

We present a detailed theoretical analysis of very rare, exclusive hadronic decays 
of the electroweak gauge bosons V = W, Z from first principles of QCD. Our main 
focus is on the radiative decays V —My, in which M is a pseudoscalar or vector 
meson. At leading order in an expansion in powers of A^cD/my the decay amplitudes 
can be factorized into convolutions of calculable hard-scattering coefficients with the 
leading-twist light-cone distribution amplitude of the meson M. Power corrections to the 
decay rates arise first at order [Aqcd/ my)'^. They can be estimated in terms of higher- 
twist distribution amplitudes and are predicted to be tiny. We include one-loop 0{as} 
radiative corrections to the hard-scattering coefficients and perform the resummation 
of large logarithms (os ln(my//iQ))"' (with po ~ 1 GeV a typical hadronic scale) to 
all orders in perturbation theory. Evolution effects have an important impact both 
numerically and conceptually, since they reduce the sensitivity to poorly determined 
hadronic parameters. We present detailed numerical predictions and error estimates, 
which can serve as benchmarks for future precision measurements. We also present an 
exploratory study of the weak radiative decays Z —)• MW. Some of the decay modes 
studied here have branching ratios large enough to be accessible in the high-luminosity 
run of the LHC. Many of them can be measured with high accuracy at a future lepton 
collider. This will provide stringent tests of the QCD factorization formalism and enable 
novel searches for new physics. 


1 Introduction 


One of the main challenges to particle physics is to obtain a rigorous control of strong- 
interaction phenomena in a regime where QCD is strongly coupled. Over the years, lattice 
QCD has made much progress in computing the static properties of hadrons from hrst princi¬ 
ples. The concept of quark-hadron duality has enabled us to make systematic predictions for 
inclusive decay processes with a large energy release, such as hadrons at large or 

inclusive weak decays like B Xlu. In these cases, non-perturbative aspects of the strong 
interactions can be accounted for using a local operator-product expansion. A conceptually 
more difficult problem is to control strong-interaction effects in exclusive hadronic processes at 
large energy. For deep-inelastic scattering a factorization theorem can be derived, in which all 
non-perturbative physics associated with the initial-state nucleon can be described in terms of 
parton distribution functions (PDFs), up to power corrections suppressed by Aqcd/\/s- The 
same framework is routinely used to calculate cross sections at hadron colliders such as the 
LHC in terms of convolutions of calculable partonic cross sections with PDFs, even though 
the underlying factorization formula can only be proved for the simplest such processes. 

The QCD factorization approach developed by Brodsky and Lepage 0,3 , Efremov and 
Radyushkin 3 3 and others 3 provides a theoretical basis for controlling strong-interaction 
effects in exclusive processes with individual, highly energetic hadrons in the final state. 
Bound-state effects are accounted for in terms of light-cone distribution amplitudes (LCDAs) 
of these hadrons, which are defined in terms of the matrix elements of non-local quark and 
gluon operators with light-like separation. This approach provides an expansion of amplitudes 
in powers of Aqcd/Q, where Q is the large energy released to the hadronic final state. While 
the leading term can be calculated in a model-independent way, it is generally not guaranteed 
that power corrections can be meaningfully computed, as they may involve ill-defined overlap 
integrals. About 15 years ago the QCD factorization formalism was generalized to deal with a 
particularly complicated class of processes: non-leptonic, exclusive weak decays of B mesons 
[6|-l9|. The additional complication consists in the presence of soft form-factor contributions 
in addition to the hard-scattering contributions described in the Brodsky-Lepage framework. 
The decay amplitudes are expanded in powers of Aqcd/E, where E nib is the energy of the 
final-state hadrons in the R-meson rest frame. While the predictions obtained at leading order 
are theoretically clean and in reasonable agreement with experiment, already the hrst-order 
power corrections involve ill-dehned overlap integrals. This introduces poorly known model 
parameters and makes phenomen olog ical predictions less precise. Since the advance of soft- 
collinear effective theory (SCET) |l0l-ll3l|. the QCD factorization approach can be rephrased 
in the language of effective held theory, which helps making its workings more transparent. 

All existing applications of the QCD factorization approach have suffered from the fact 
that the characteristic energy scales are not sufficiently large for power corrections to be 
negligible (see e.g. 0 for a recent discussion). It is then notoriously difficult to disentangle 
power-suppressed effects from the uncertainties related to the shapes of the hadron LCDAs. 
Unfortunately, no comprehensive experimental program to determine the leading-twist LCDAs 
of the ground-state mesons and baryons - analogous to the large-scale effort to determine the 
PDFs of the proton with high accuracy - is conceivable. In this paper we propose using 
exclusive decays of the heavy electroweak gauge bosons W and Z into final states containing 
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a single meson as a laboratory to test and study the QCD factorization approach in a context 
where power corrections are definitely under control. The enormous rates of W and Z bosons 
that will become available at future colliders will present us with a new playground for precision 
electroweak and QCD physics, which will make such studies feasible. With 3000 fb“^ collected 
during the high-luminosity run at the LHC, one will have produced more than 10^^ Z bosons 
and 5 • 10^^ W bosons in both ATLAS and CMS. A clean, tagged sample of W bosons is 


expected to come from top-quark decays [l5(]. At a future lepton collider such as TLEP, 


samples of up to 10^^ Z bosons per year can be expected in a dedicated run at the Z pole 
(l^ . Our main focus is on the simplest processes: the hadronic radiative decays Z —> My and 
W —)■ My, where M is a pseudoscalar or vector meson. They offer a perfect way to probe some 
properties of the leading-order LCDAs of various mesons. The price one needs to pay is that 
the higher the energy release in the process, the smaller the probabilities for any particular 
exclusive final state are. The branching fractions we obtain range from few times 10“® to few 
times 10“^^ or even smaller. The big challenge of such a program will be to measure such 
decays experimentally with some precision. While we do not perform a detailed feasibility 
study in this work, we speculate that some of these rare modes will be accessible in the high- 
luminosity run of the LHC, at a level that will be useful to probe our theoretical predictions. 
At future lepton colliders operating on the Z pole, it would be possible to measure several 
of these decays at or below the 1% level. This may present us with a unique opportunity to 
extract information about LCDAs in a theoretically clean environment. 

Our interest in this subject was raised by recent investigations of the exclusive decays 
h ^ V'y 17-20] and h —)■ ZV j^, ^ of the Higgs boson to hnal states containing a single 


vector meson. It was proposed to use these decays as a way to probe for possible non-standard 
Yukawa couplings of the Higgs boson. Such measurements are extremely challenging at the 
LHC and other future colliders. Observing exclusive hadronic decays of W and Z bosons would 
provide a proof of principle that this kind of searches can be performed. An encouraging first 


search for the decays Z —)■ J/V'q and Z —)■ T(nS') y has just been reported by ATLAS [23 


From a theoretical perspective, the very rare, exclusive radiative decays of W and Z 
bosons have received relatively little attention in the literature, and very few accurate pre¬ 
dictions for such branching fractions have been obtained. In a pioneering study jl^, Arnelos, 
Marciano and Parsa presented a hrst detailed analysis of the decays W —)■ Py and Z —)■ Py for 
both light and heavy pseudoscalar mesons in the final state. Strong-interaction effects were 
parametrized in terms of vector and axial-vector form factors, which were estimated using 
ideas from perturbative QCD on the asymptotic behavior of form factors at large momentum 
transfer. Several years later, Manohar studied the decays Z —)■ irW and Z ^ using a 


local operator-product expansion [2^, which expresses the decay amplitudes as power series 
in parameters Uq = 0.26 and Uq = 2, respectively. If only the leading term is 

kept, the amplitudes can be related to the pion matrix element of the axial-vector current, 
which is proportional to the pion decay constant /^r. For radiative decays such as Z —)■ vry this 
truncation cannot be justified theoretically, and the infinite tower of local operators would 
need to be resummed. In a very recent work, the method developed by Manohar was used to 
derive an estimate for the IF —?• vry branching fraction [l^. While the analyses performed in 
these papers can provide some order-of-magnitude results, they do not allow to obtain accu- 
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rate predictions with reliable error estimates. In a classic paper |2^ , Guberina et al. analyzed 
the radiative decays of the Z boson into heavy qnarkonia in the non-relativistic limit. The 
first relativistic corrections to the Z —)■ and Z —)■ T(1S')7 decay rates were added only 

recently in j^, where in addition the anthors considered for the first time the decay Z —>■ 07, 
using an approach closely related to ours. As we will discuss, renormalization effects have a 
profound impact on the decay amplitudes. When evolved up to the relevant scales of order 
the Z-boson mass, the LCDAs of heavy qnarkonia can no longer be accurately described by 
the leading term in a non-relativistic expansion. 

In the present work, we present a comprehensive analysis of a large class of radiative 
decays of W and Z bosons using the QCD factorization approach, including for the hrst 
time a consistent treatment of 0{as) corrections and performing the resummation of lar^ 
logarithms of order ln(m|//ig)) , with /ig ~ 1 GeV, to all orders in perturbation theory^ 
Our approach provides a systematic expansion of the decay amplitudes in powers of the small 
parameters asijnz) ~ 0.1 and Agcn/’^z ~ 0.01. We study the structure of the leading power 
corrections to the Z —)• My and W —)■ My decay rates and show that they are of second 
order and hence negligibly small, of order 10“^ relative to the leading terms. For processes 
involving heavy quarks, power corrections of order {mq/mzY exist, which are still very small 
(less than 1%) even for hnal-state mesons containing 6-quarks. Finally, using the most recent 
experimental data we perform a reanalysis of meson decay constants, which provide crucial 
input to our phenomenological analysis. 

Our paper is organized as follows: In Section [2] we derive a factorization theorem for the 
Z —)■ My and W —)■ My decay amplitudes, in which they are expressed as convolutions of 
calculable hard-scattering kernels with meson LGDAs. We explain how the kernel functions 
can be calculated by performing projections of on-shell partonic amplitudes. We then summa¬ 
rize the existing theoretical information on the shapes of the LGDAs for both light and heavy 
mesons and study their behavior under scale evolution. In Section [3] we apply this approach to 
derive explicit predictions for the Z —)■ My and W —)■ My decay amplitudes at leading power 
in Aqcd/ mz- The relevant convolution integrals of hard-scattering kernels with LGDAs are 
calculated in analytic form using an expansion in Gegenbauer polynomials. We demonstrate 
that renormalization-group (RG) evolution from a low hadronic scale up to the electroweak 
scale of relevance to these processes has the nice effect of signihcantly reducing the sensitivity 
to poorly determined hadronic parameters. By studying radiative decays into transversely 
polarized vector mesons, we present some detailed estimates of power-suppressed effects in 
the QGD factorization approach. In some old papers, it was suggested that the radiative 
decay amplitudes into pseudoscalar mesons can be hugely enhanced due to effects of the axial 


anomaly j28l. l29[| . We explain why such an enhancement does not exist. We then present our 


numerical predictions for VF, Z —)■ My branching fractions, including detailed error estimates. 
The results span more than three orders of magnitude, and we explain the striking differences 
seen between the various decay channels in terms of electroweak couplings, differences in decay 
constants, and enhancement factors occurring for heavy-light mesons. In Section 0] we present 
a hrst exploratory study of the weak radiative decays Z —)■ MW, in which a heavy W boson is 


^Radiative corrections to the Z —>• J/tp^ and Z —>• T(1S')7 decay amplitudes were included in [27| in the 
non-relativistic limit, but they were not included so far in any analysis of decays into light final-state mesons. 
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Figure 1: Leading-order Feynman diagrams for the radiative decays The 

meson bound state is represented by the gray blob. 


part of the hnal state. In this case signihcantly less energy is released to the hnal-state meson 
M, and as a result the QCD factorization approach can be tested at energies of order 10 GeV, 
about a factor 2 higher than those available in exclusive S-meson decays. We round off our 
study in Section O with some experimental considerations. Our main results are summarized 
in Section [6l Technical details of our calculations and the extraction of meson decay constants 
are relegated to three appendices. 


2 Theoretical framework 

Our main focus in this work is on the rare, exclusive radiative decays Z —)■ My and W —)■ My, 
where M denotes a pseudoscalar or vector meson. We assign momentum k to the hnal-state 
meson and q to the photon. The leading-order Feynman diagrams for the case of Z —)■ My 
are shown in Figure [H The decay plane is spanned by the vectors k and q. We will refer to 
vectors in this plane as being longitudinal, and to vectors orthogonal to it as being transverse. 
We only consider cases where the mass of the hnal-state meson satishes tum ^ fnz- Up to 
corrections suppressed as {mM/^zYi this mass can then be set to zero. In this limit, we have 
k^ = En^ and q^ = En^, where E = mz/^ is the energy of the hnal-state particles in the 
Z-boson rest frame, and n and n are two light-like vectors satisfying n ■ n = 2. 


2.1 Derivation of the factorization formula 


For the purposes of this discussion we work in the rest frame of the decaying heavy boson. The 
decay amplitudes can be calculated from hrst principles using the QCD factorization approach 
because the energy E released to the hnal-state meson is much larger than the scale 
of long-distance hadronic physics. At leading power in an expansion in Aqcd/'^z, they can 
be written as convolutions of calculable hard-scattering coefficients with LCDAs of the meson 
M. A s imp le way to derive the corresponding factorization theorem employs the formalism of 


SCET |inl-ll3|. It provides a systematic expansion of decay amplitudes in powers of a small 


expansion parameter A = Aqcd/E. The light hnal-state meson moving along the direction 
can be described in terms of collinear quark, anti-quark and gluon helds. These particles carry 
collinear momenta Pc that are approximately aligned with the direction n. Their components 
scale like (n -pc, h -Pc^pj:) ~ E(A^, 1, A). Note that p^ ~ Aq(-.j 3 , as appropriate for an exclusive 
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hadronic state. The collinear quark and gluon fields are introduced as gauge-invariant objects 
dressed with Wilson lines. Explicitly, one defines 


30, 31 




where iD^ = id^ + denotes the covariant collinear derivative, and 

r-O 


Wc{x) = Pexp j 


dtn ■ Ac{x + tn) 


( 1 ) 

( 2 ) 


is a collinear Wilson line extending from x to infinity along the direction n. Both fields are 
of 0{\) in SCET power counting. Adding more component fields to an operator always leads 
to further power suppression. At leading order in A, the operators with a non-zero matrix 
element between the vacuum and a single meson state are thus of the form Xc{tn) ... A’c(O) 
and A'^j_(tn) ... Alc±/i(0), where without loss of generality we set a: = 0 for one of the fields. 
Since the effective collinear fields are gauge invariant by themselves, composite operators built 
out of these fields can be non-local along the light-like direction n. The two-gluon operator 
would only be relevant for decays into mesons containing a flavor-singlet component on their 
wave functions, such as the pseudoscalar mesons rj and rj' [^. Such decays will be discussed 
in a forthcoming publication js^. It follows that at leading power in the expansion in A, the 
Z —)■ M'j and W —)■ My decay amplitudes into non-singlet final states can be written in the 
factorized form 



— t/i 

dtCi{t, fx) {M{k) \ Xcitn) — Tj A’c(O)lO) -|- power corrections 
dt Ci{t, fi) {M{k) \ q{tn) — Tj [tn, 0] q'(0)|0) -|- power corrections. 


( 3 ) 


where fx is the factorization scale, and Tj G { 1 , 75 , 7 ^}. The four matrices (^/2) Tj provide a 
basis of Dirac matrices sandwiched between two collinear quark spinors. The Wilson coeffi¬ 
cients Ci(t) are process dependent and can be calculated perturbatively. In the last step we 
have used the definition ([1]) and combined the two Wilson lines Wdtn) W^{0) = [th, 0] into a 
straight Wilson line extending from 0 to tn. The meson matrix elements of the bi-local oper¬ 
ators in the second line define the leading-order LCDAs of pseudoscalar and vector mesons. 
Specifically, one has 


{M{k) I q{tn) | ( 75 ) [tn, 0] g(0) |0) = -ifnE dx /i); M = P, Vj|, 

(I4(A;)| g(M) ^ 7 ^ [fh, 0] g(0)|0) = da; 0^(a;,/i), 


( 4 ) 


where E = n ■ k/2 denotes the energy of the meson in the rest frame of the decaying boson, 
fp and fv are the decay constants of pseudoscalar and vector mesons defined in terms of their 
matrix elements of local (axial-)vector currents, and fy{g) is a scale-dependent vector-meson 
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decay constant defined in terms of a matrix element of the QCD tensor current. The leading- 
order LCDAs can be interpreted as the amplitudes for finding a quark with longitudinal 
momentum fraction x insinde the meson. The factor of 75 in the first equation is present for a 
pseudoscalar meson (M = P) but absent for a longitudinally polarized vector meson (M = V||). 
The projection onto a transversely polarized vector meson does not arise at leading power in 
the radiative decays of W and Z bosons. For a given meson, exactly one of the possible 
Dirac structures contributes, and we denote the corresponding Wilson coefficient by 
Dehning the Fourier-transformed Wilson coefficient, called the hard function, via 

Hm{x, = Jdt CM{t, fi) , (5) 


we obtain the factorization formula 


A = —ifME / dx Hm{x, /i) 4>m{x, /i) + power corrections . 


( 6 ) 


Insertions of additional collinear fields or derivatives yield power-suppressed contributions. In 
particular, the insertion of an additional collinear gluon field gives rise to three-particle LCDAs. 
In order to fully establish the factorization theorem ([2]) one must show that the convolution 
integral over the momentum fraction x converges at the endpoints. This question has been 


addressed in the context of the more complicated processes B —)■ 7 /z/ [3^ and B —)■ iF *7 


m 


The behavior near the endpoints can be described by means of soft-collinear fields 
3^,1^ with momenta scaling as (n-pscn-pscPic) ~ F^(A^, A, The contributions of such 


modes are always power suppressed. In the present case, we hnd that endpoint singularities 
are absent at leading and subleading power in the large-energy expansion. 

LCDAs play the same role for hard exclusive processes which PDFs play for inclusive 
ones. While they encode genuinely non-perturbative hadronic physics, they can be rigorously 
defined in terms of non-local operator matrix elements in QCD [iH^. These matrix elements 
can be systematically expanded in terms of structures of different twist. When applied to high- 
energetic exclusive processes such as the ones considered here, the twist expansion translates 
into an expansion in powers of Aqqd/E. There is an extensive amount of literature devoted 
to the study of distribution amplitudes. For light pseudoscalar mesons, the two- and three- 
particle LCDAs up to twist-3 order were studied, e.g., in [^, while the corresponding LCDAs 


for vector mesons were analyzed, e.g., in |39l-l4l[|. We stress that, at the scale of the large 


energies released in decays of W and Z bosons, even charm and bottom quarks can be treated 
as light quarks, and hence heavy mesons containing these quarks can be described by LCDAs. 
This will be discussed further below. 

In order to apply these results in practical calculations, it is convenient to define momentum- 
space projection operators, which can be applied directly to the decay amplitudes computed 
with on-shell external parton states j^, ^]. For all two-particle projections onto LCDAs of 
leading and subleading twist, it is sufficient to assign momenta ki = xk + k± + ... and 
k 2 = {1 — x)k — k± + ... to the quark and the anti-quark in the meson M, where k is treated 
as a light-like vector (/c^ = 0). Meson mass effects of order enter only at twist-4 level. 
They have a tiny numerical impact for the decays considered here, and we will consistently 
set —)■ 0 unless noted otherwise. The variables x and (1 — x) denote the longitudinal 
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momentum fractions carried by the quark and the anti-quark in the two-body Fock state of 
the meson. Each Feynman diagram gives an expression of the form 


u{ki) A{q, fci, fca) v{k 2 ) = Tr [^(^ 3 ) ^(^ 1 ) A{q, fci, ^ 2 )] 


-)■ 


dxTi [MM{k,x,^i) A{q, ki, k 2 )\ 


(7) 


>0 > 


where in the last step we have introduced the light-cone projection operator MM{k,x, fi) for 
the meson M, which at higher order contains derivatives with respect to the parton transverse 
momentum k±. It is understood that k± is set to zero after these derivatives have been 
performed. 

Up to twist-3 order, the light-cone projector for a pseudoscalar meson can be written in 


the form [ 8 |, ^2 


^ f 

Mp{k, x,fx) = ^75 0 p(a;, p) - /ip(/i) 75 


(j)p(^Xj fji') icTp 


k^rA (f)'^{x,^) 


+ iaf,^k^ 




6 


dk 


±u 


k ■ n 6 
-|- 3-particle LCDAs > . 


( 8 ) 


Here 0p is the leading-twist LCDA of the meson, while (j)p and (j)a- denote the two-particle 
LCDAs appearing at twist-3 order. These are scale-dependent functions, which we define in 
the MS renormalization scheme. The decay constant fp of the meson P is defined in terms of 
its matrix element of a local axial-vector current 


{P(/c)| |0) = -ifpk^ . 


(9) 


The scale-dependent parameter governs the normalization of 

the twist-3 LCDAsH The vector n in the above expression denotes a longitudinal light-like 
vector not aligned with k. A convenient choice is to take the photon momentum, n = q. 
At twist-3 order the projector also contains three-particle LCDAs containing a quark, an 
anti-quark and a gluon. We will see that the contributions of twist-3 LCDAs are strongly 
suppressed compared with those of the leading-twist amplitudes. In order to estimate their 
effects, we will for simplicity neglect the three-particle LCDAs. This is referred to as the 
Wandzura-Wilczek approximation (WWA) jd^. When this is done, the QCD equations of 
motion fix the form of the twist-3 LCDAs completely, and one obtains js^ 

(f)p(^x, p.) 1 , p) 6 x(l x). ( 10 ) 


The light-cone projection operators for vector mesons are more complicated. They are given in 
Appendix IS For our purposes it suffices to quote the projector for a longitudinally polarized 
vector meson at leading power. It is 


(/C, X, p) = ^ (j)y(^x, p) + • • • = ^ (x, p) + . . . 

" A k • n 4 


^Note that holds for charged and neutral pions, see e.g. [32 1. 


( 11 ) 
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The function 0y(a:,/i) is sometimes called in the literature. We have used that the 

longitudinal polarization vector is given by — 'my f^)- The vector-meson decay 

constant /y is dehned in terms of the local matrix element 

(V(k,ev)lgi7''Q2 |0) = -ifvmve*v ■ (12) 

Before proceeding, let us comment on the strnctnre of power corrections to the factorization 
formula (jS]). Inspecting the explicit form of the projection operator for a psendoscalar meson 
in (jH]), and the corresponding projectors for vector mesons given in flA.lD and flA.4D . we observe 
that consecutive terms in the twist expansion contain even and odd nnmbers of Dirac matrices 
in alternating order. Since the gange interactions in the Standard Model preserve chirality, 
it follows that for a given helicity amplitude either all terms with an even nnmber of Dirac 
matrices contribnte or all terms containing an odd nnmber, bnt not both. Conseqnently, the 
SCET expansion for the Z —)■ My decay amplitndes with hxed polarizations of all particles is 
an expansion in powers of (Aqcd/'^z)^- The power connting changes when qnark-mass effects 
are taken into acconnt. They give rise to chirality-changing vertices, which give corrections 
snppressed by mg/m^ to both the amplitndes and the meson projectors. This leads to power 
corrections of order mgAgco/'^l and {mq/mzY■ For heavy qnarks with mg S> Aqcd, the 
latter corrections are the dominant ones. However, as long as the relevant qnark masses mq 
are much smaller than the hard scale of the process, these corrections are still small. The 
present case is different from the sitnation encountered in exclusive H-meson decays [^-[^, 
where the presence of a heavy qnark mass, which is of the same order as the energy released in 
the decay, allows for 0{1) chirality-changing interactions. In this case the decay amplitudes 
receive hrst-order KqQjy/rrib corrections. 


2.2 Systematics of the Gegenbauer expansion 

The leading-twist LCDAs obey an expansion in Gegenbaner polynomials of the form 0,3 


4>m{x, /i) = 6a:(l — x) 




71=1 


(13) 


which can be inverted to give 




(h) = 


2(2n + 3) 
3(n + l)(n + 2) 


dxC^^‘^\2x - 1 ) 4>m{x,iZ) . 


(14) 


The Gegenbauer moments have a diagonal scale evolntion at leading order in pertnrbation 
theory. They are non-perturbative hadronic parameters, which can only be accessed using 
data or a noii-pertnrbative approach snch as light-cone QGD snm rules (see e.g. j^-M|) or 
lattice QGD . In Table [U we collect the values for the decay constants and the hrst two 
Gegenbaner moments a ^2 light pseudoscalar and vector mesons. Our notation is such that 
K(*) (qs) with q = u,d, and x is the momentnm fraction of the light qnark q. 

An expansion such as ffT5D is useful provided we have some reason to believe that the 
inhnite series is dominated by the hrst few terms. Higher-order Gegenbaner moments of the 











Meson M 

fM [MeV] 

af (/io) 

a^ipo) 

TT 

130.4 ±0.2 

0 

0.29 ±0.08 

K 

156.2 ±0.7 

-0.07 ±0.04 

0.24 ±0.08 

P 

212 ±4 

0 

0.17 ±0.07 

U! 

185 ±5 

0 

0.15 ±0.12 

K* 

203 ±6 

-0.06 ±0.04 

0.16 ±0.09 

0 

231 ±5 

0 

0.23 ±0.08 


Table 1: Hadronic input parameters for light pseudoscalar and vector mesons, with 
scale-dependent quantities dehned at /xq = 1 GeV. We assume isospin symmetry and 
use the same values for charged and neutral mesons. The values for /^r and Jk are taken 
from I45II . The other decay constants are extracted from t~ — )■ M~Uj- and — )■ 

decays , as discussed in Appendix |Bl For all other parameters we adopt the values 
compiled in from a combination of results obtained using lattice QCD ji^ and 
light-cone QCD sum rules (see e.g. (iol. Idsl-I^b including conservative error estimates. 


pion were studied in [52, l53| using a QCD sum-rule approach employing non-local vacuum 


condensates. These authors hnd = 0.20, aj = —0.14, Og = 5 • 10“^, and Og = aJo = 4 • 10“^ 
at the scale fio = 1 GeV. Their value of aj is consistent with the result given in Table [H while 
higher moments with n > 6 are estimated to be negligibly small. On the other hand, in 
more recent work ^ the authors have performed hts to the hrst eight Gegenbauer moments 
of the pion LCDA using data on the 7r°7*7 form factor obtained by the BaBar and Belle 


collaborations [55j, [5^. They hnd aj = 0.10(0.14), aj = 0.10(0.23), Og = 0.10(0.18) and 


Og = 0.034(0.050) at juq = 1 GeV for Belle (BaBar), which suggests that Og and Og may not 
be insignihcant. In our phenomenological analysis we will vary a^(/io) between —0.15 and 
-1-0.15 for all light mesons and use this to estimate the effect of unknown higher Gegenbauer 
moments. With this treatment, the relevant combination of Gegenbauer coefficients given in 
relation fl44jl below agrees with all of the above models within our quoted uncertainties. 

It is an important question to ask what can be said on general grounds about the behavior 
of the Gegenbauer expansion. It is commonly assumed, and is supported by power-counting 
analyses in SCET, that the leading-twist LCDAs vanish at the endpoints x = 0 and x = 1, 
such that the integrals ^ </>m(x) and converge. This statement implies that 

the inhnite sums Yhn^n converge. Barring accidental cancellations, this 

requires that for large n the coefficients fall off faster than l/n, and this condition should 
hold for all values of uq- From a physical point of view, high-rank Gegenbauer polynomials 

/o /o\ 

Cn ' T2x — 1) with n S> 1 resolve structures on scales Ax ~ 1/n. For a light meson M, it is 
reasonable to assume that the LGDA 0 m ( 2 ^) does not exhibit pronounced structures at scales 
much smaller than 0{1), in which case the coefficients must decrease rapidly at large n. 

The LGDAs of heavy mesons are an exception to this rule, since the presence of the heavy- 
quark mass introduces a distinct scale. For a quarkonium state M ~ {QQ) composed of two 
identical heavy quarks, the LGDA peaks at x = 1/2 and has a width that tends to zero in the 
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limit of infinite heavy-quark mass. The second moment of the LCDA around x = 1/2 can be 



1 / 2 \ 
dx (2x - 1)^ /io) = —+ 0{v‘^). 


(15) 


To derive this result one uses that in the heavy-quark limit x = 


PQ-^ 

2moV-n 


— 1+Vz 
2 


where is 


a light-like vector, and pq = mqV^ + denotes the momentum of the heavy quark inside 
the quarkonium state with velocity V^. The various vectors are defined such that V ■ n = 1 
and V ■ k = D. In the rest frame of the quarkonium state we can choose = (1,0), 
= (1, —e^), and k^ = {0,mQv), where the 3-vector v is the residual velocity of the heavy 
quark inside the (QQ) bound state. The factor 1/3 on the right-hand side of flTHll is due 
to rotational invariance in the rest frame. Numerical values for the NRQCD matrix element 
(n^) for the J/ip and T(IS') states have been obtained from an analysis of the leptonic decay 
rates —)■ e'''e“) and r(T(lS') —)■ e'''e“) including hrst-order corrections and non- 

perturbative contributions proportional to In this way, the values (n^)j/^ = 0.225/g;Qgg 
jool l and {v‘^)t(is) = —0.009 ±0.003 have been extracted, the latter one being inconsistent 
with the fact that the second moment in f[T5|l must be positive. Both estimates suffer from 
the fact that the two-loop j^, 631 and three-loop perturbative corrections to the NRQCD 
predictions for these decay rates are known to be huge, precluding a reliable extraction of non- 
perturbative parameters. Based on the power-counting rules of NRQCD one would naively 
expect that (n^)j/^ ~ 0.3 and (n^)T(is) ~ 0.1, and we will use these estimates, along with a 
50% relative error assigned to them, in our phenomenological analysis. For our calculations 
we need the hrst inverse moments of the LCDA with respect to x or (1 — x). Expanding the 
inverse moments about x = 1/2, it is immediate to derive the model-independent relation 


4>m{x,po) ^ (j)M{x,po) ^ 2 

Io X Jo 1 - a; 


1 + + o(.„4) 


( 16 ) 


As a reasonable model at the low scale po = 1 GeV we adopt the Gaussian ansatz 


(pnix, iio) = N„ exp 


\/^i 


TTCT 


X 


1\2 


2 cr 2 


2 _ \v‘^)m 


a = 


12 


(17) 


where the polynomial in front of the Gaussian factor ensures that the LGDA vanishes at the 
endpoints x = 0,1. The normalization constant ^ 1 can be expressed in closed form in 
terms of an error function. 

For a heavy-light meson state M ~ {qQ) composed of a light quark and a heavy anti-quark, 
the LGDA peaks at a small value x ~ Aqcd/’^m, where x refers to the momentum fraction of 
the light spectator quark. The appropriate effective field theory for heavy-light bound states is 
called heavy-quark effective theory (HQET), see for a review. In the context of this theory, 
it is possible to show that the hrst moment of the LGDA is determined by the ratio Am/ptim, 
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Meson M 

/m [MeV] 

Am [MeV] 


a 

D 

204.6 ±5.0 

460 ± 110 

- 

0.246 ± 0.059 

Dg 

257.5 ±4.6 

550 ± 150 

- 

0.279 ± 0.076 

B 

186 ±9 

460 ± 110 

- 

0.087 ±0.021 

Bs 

224 ± 10 

550 ± 150 

- 

0.102 ±0.028 

JlJj 

403 ±5 

- 

0.30 ±0.15 

0.158 ±0.040 

T(15) 

684 ±5 

- 

0.10 ±0.05 

0.091 ±0.023 

T(45) 

326 ± 17 

- 

0.10 ±0.05 

0.091 ±0.023 


Table 2: Hadronic input parameters for pseudoscalar and vector mesons containing 
heavy quarks. Scale-dependent quantities are defined at /xq = 1 GeV. The values for 
fn and f^s are taken from ji^. The values for fs and /b are taken from two recent, 
unquenched lattice calculations 00, which obtain identical central values but quote 
very different error estimates. We quote the averages of the uncertainties given by the 
two groups. The values of the J/ijj and T(nS') decay constants can be derived from 
data, as explained in Appendix iBl 


where tum denotes the heavy-meson mass and Am = mM — mg (with mg being the pole 
mass of the heavy quark) is a hadronic parameter. One obtains (x) = ^Am/thm + 0[as{mg)] 
0 ], where the one-loop radiative corrections have been calculated in 0] and are numerically 
significant. In our analysis below we need the first inverse moment of the LCDA with respect 
to X, which is of order ruM/AgcD and cannot be related to a local HQET matrix element. 
One defines 

Jo ^ Am(/^o) 

where the hadronic parameter Am(/^o) ~ Aqcd is independent of the heavy-quark mass, and 
the dots denote corrections that are power-suppressed relative to the leading term. The 
parameter Am is poorly known at present. A QCD sum-rule estimate for the B meson yields 
As(l GeV) = (460±110) MeV [0], and we will use this value in our phenomenological analysis 
for both B and D mesons. Goncerning Bg and Dg mesons, we shall use the estimate Am^ ~Am ~ 
90 MeV from 0] and increase the error to ±150 MeV. As a plausible model at a low scale 


/io = 1 GeV we take 




x(l — x) 


a 


2 


exp 



Am(/^o) 

mM 


(19) 


where the normalization constant ^ 1 can be determined in closed form. For heavy- 
light mesons M ~ {Qq) containing a heavy quark and a light anti-quark, one simply replaces 
X -f-)- (1 — x) in the above relations. 

In Table [2] we collect the values for the decay constants and the width parameters for 
heavy pseudoscalar and vector mesons, which will be used in our phenomenological analysis. 
In the cases of {qQ) and {QQ) bound states, Gegenbauer moments of roughly n < 1/a 
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give important contributions to the LCDAs, because they are required to resolve the narrow 
structures of the LCDAs near the peak region. For example, at po = 1 GeV the hrst 5 ( 6 ) 
Gegenbauer coefficients of the i?-meson (T-meson) LCDA are larger in magnitude than 0.1, 
and the hrst 7 (12) Gegenbauer coefficients are larger than 0.01. We will discuss in the next 
section that the effects of QGD evolution from a low scale /ig up to a high scale reduces the 
high-rank Gegenbauer moments much stronger than Gegenbauer moments of low rank. For 
example, at /i = mz only the hrst 3 (2) Gegenbauer coefficients of the i?-meson (T-meson) 
LGDA are larger than 0.1, and the hrst 6 ( 8 ) Gegenbauer coefficients are larger than 0.01. 
As a result, the shapes of the LGDAs for mesons containing heavy quarks are signihca ntly 
ahected by RG evolution. For the case of the R-meson LGDA, this ehect was studied in [72|. 
Gonsequently, the low-scale predictions for the inverse moments considered here are strongly 
modihed at /x = 0{mz)- 


2.3 Radiative corrections and RG evolution 


In order to improve the accuracy of our predictions and be in a position to meaningfully 
discuss the setting of the factorization scale /x, we include the 0{as) radiative corrections 
to the leading-twist contributions in our analysis, hnding that RG evolution ehects are very 
important. The reason is that logarithms of the form (os ln(m|//iQ)) , where /xq ~ 1 GeV 
denotes the scale at which non-perturbative calculations of the LGDAs are performed, are 
numerically large and must be resummed to all orders of perturbation theory. We perform 
the calculation of the loop diagrams shown in Figure fusing dimensional regularization with 
d = 4 — 2e space-time dimensions. The individual on-shell graphs contain both UV and IR 
divergences. For the decays Z —)■ My and W —)■ My, which are mediated by vector and 
axial-vector currents, the UV divergences cancel in the sum of all diagrams. The remaining 
IR poles are cancelled when we renormalize the LGDAs. To this end, we express the bare 
LGDAs in terms of the renormalized ones, 

( 20 ) 

Jo 

At one-loop order, one obtains 13,0 

Z^{x, y, y.) = 5{x - y) + Vo{x, y) + 0{a ^,), (21) 

where Cp = {N^ - l)/{2Nc) = 4/3, and 


Vo{x,y) 


Ad-y) 

x{l - y) 

X , , 1 — X 

- 

e{y x) + 

d(x - y) 

[.y 1 - 1 / 



- x) 
y-x 


+ yii 


x) 


0{x - y) 
x-y 


( 22 ) 


is the one-loop Brodsky-Lepage kernel. For symmetric functions g{x,y), the plus distribution 
is dehned to act on test functions /(x) as 


dy[g{x,y)] f{x)= / dy g{x,y) [f{x) - f{y)] . 


(23) 
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Figure 2: One-loop QCD corrections to the first diagram in Figure [U Analogous 
corrections exist for the other diagram. 


Besides the subtraction of 1/e poles using dimensional regularization in the MS scheme, 
one must carefully address the question of how to define 75 in d 7 ^ 4 dimensions. Some of 
the amplitudes considered in this work involve traces of Dirac matrices containing a single 
insertion of 75. It is well known that for such traces the naive dimensional regularization 
scheme with anti-commuting 75 is algebraically inconsistent. Here we employ the ’tHooft- 
Veltman (HV) scheme 73(], in which 75 = i 7 ° 7 ^ 7 ^ 7 ^ anti-commutes with the four matrices 7 ^ 
with p G {0,1, 2, 3}, while it commutes with the remaining (d — 4) Dirac matrices 7^11 While 
this dehnition is mathematically consistent, it violates the Ward identities of chiral gauge 
theories by hnite terms, which must be restored order by order in perturbation theory [7^. In 
the present case, this is accomplished by performing the hnite renormalization = ^hv^hv 
of the axial-vector current, where 


Z„v(fi) = 1 - AAd + 0(7). (24) 

71 

In addition, the leading-twist LCDA of a pseudoscalar meson, which is dehned in terms of a 
matrix element of a non-local axial-vector current on the light-cone, receives a hnite renor¬ 
malization of the form 


where 


77 


h) — 


dyZ^y{x,y,fi) (j)p{y,y) , 


(25) 


ZuU^,y,y) = d{x -y) + 


2C'i?Q!s(p) 


TT 


- d{y -x) + - - 9 {x - y) 

y 1 -d 


+ • (26) 


This redehnition is important to restore the proper normalization of the LCDA 0p(a:,p). 
Integrating relation fl25|l over x, we hnd that 


/ dx0p,Hv(a;,/i) = ^Hv(/^) / dy(j)p{y,n) = ^hv(/^) > 
10 Jo 


(27) 


with Zhv given in fl24j) . The integral turns the matrix element of the non-local axial-vector 
current into the corresponding local matrix element. 


^For the purposes of our analysis, 
proposed by Larin [74 1. 


the HV scheme is equivalent but more convenient than the scheme 
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Our final expressions for the decay amplitudes will contain the scale-dependent, leading- 
twist LCDAs 0A4-(x,/i) with M = P,V\\. These functions satisfy the integro-differential evolu¬ 
tion equation 


h -r- 

djjL 


dyV{x,y,fx) (j)M{y,y), 


(28) 


where V{x,y,ix) = Vo(x, y) _j_ C>(q/ 2)_ xhe eigenfunctions of the one-loop Brodsky- 

Lepage kernel Vo{x,y) in are the Gegenbauer polynomials 6a:(l — x) Cn^^\2x — 1), and 
hence the Gegenbauer moments an{u) defined in fll3p are multiplicatively renormalized at this 
order. They obey the RG equation pj 


d M/ \ ois{y) Mt \ 

^ On (h) = -In On (h) , 


where 


In — — 


(n -h l)(n -h 2) 


-31, with if„+i = ^ 



(29) 

n+l ^ 


y-- 

^ k 

k=l 

(30) 


The evolution of the leading-twist LGDAs at two-loop order has been studied in 78l-l8l[|. The 


RG equation for the Gegenbauer moments becomes more complicated at this order, since the 
scale dependence of receives contributions proportional to with k = 0,... ,n 180- 


82| . The evolution equation can still be solved analytically using an iterative scheme. Explicit 
results for the moments up to n = 12 can be found in [1^. However, given that all present 
estimates of the hadronic parameters are afflicted with large theoretical uncertainties, it 
is sufficient for all practical purposes to use the leading-order solution fl2^ . It reads 




/ asjy) 
\ Q^s(/^o) 


7n/2^0 


(31) 


where /do = — is the hrst coefficient of the QGD /d function. Here /xq ~ 1 GeV denotes 

a low scale, at which the Gegenbauer moments are derived from a non-perturbative approach, 
while /i is a high scale to which the LGDAs are evolved. In our analysis this scale is set by the 
mass of the decaying electroweak boson. Note that one must adjust the values of /do whenever 
fi crosses a flavor threshold. All of the anomalous dimensions are strictly positive, which 
implies that —?■ 0 in the formal limit p —?• cxo. Indeed, for large n the evolution supplies 

an additional suppression factor (l/n)^ with K = In In this limit, the leading-twist 
LGDAs approach the asymptotic form 6a;(l — x). 

Figure [3] shows the RG evolution of the LGDAs of the kaon, J/'ip meson and B meson 
from a low scale po = 1 GeV up to a high scale mz- We use the Gegenbauer moments and 
width parameters collected in Tables [T] and [21 For light mesons we truncate the Gegenbauer 
expansion flT2]) at n = 2. For heavy mesons we use the model LGDAs given in fITT)) and flT2D . 
compute their first 20 Gegenbauer moments, evolve the corresponding coefficients from 
/Uo to m^, and reconstruct the LGDAs at the high scale from flT3D . The dotted line in the 
plots shows the asymptotic form 6x(l — x). Evolution effects alter the shapes of the various 
distributions in a significant way. At the electroweak scale, the LGDAs are significantly closer 
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Figure 3: RG evolution of the LCDAs of the kaon (left), the J/V’ meson (middle) and 
the B meson (right) from a low scale /iq = 1 GeV (dashed lines) to a high scale /i = mz 
(solid lines). The dotted grey line shows the asymptotic form 6x{l — x) for comparison. 


to the asymptotic form 6a;(l—a;) than at a low hadronic scale. Gonsequently, RG effects render 
our predictions more insensitive to poorly determined hadronic input parameters. Notice, in 
particular, that the LGDA of the J/ip meson at /x = mz is as close to the asymptotic form 
as the kaon LGDA. In practice, the LGDAs of heavy mesons at a scale much larger than the 
heavy-quark mass can be well described in terms of a Gegenbauer expansion truncated after 
a few Gegenbauer moments. 

2.4 Flavor wave functions of neutral mesons 

The couplings of photons and of the electroweak gauge bosons W and Z to fermions are flavor 
dependent. While the flavor content of charged mesons is unambiguous, for neutral mesons 
complications arise from the fact that a given meson can be a superposition of different flavor 
components. We write the flavor wave function of the neutral hnal-state meson in the 
form 

|V”)= E TWi ElTl"=l- (32) 

q=u,d,s,c,b q 

For heavy mesons containing charm or bottom quarks such effects can safely be neglected. 
The heavy mesons rjc and J/ip have Cc = 1, while r/h and T have Cb = 1. Mixing effects can 
however be irnportant for light mesons. 

Following |32|, we assume isospin symmetry of all hadronic matrix elements, but we dif¬ 
ferentiate between the matrix elements of mesons containing up or down quarks and those 
containing strange quarks. The 7r° and mesons are members of an isospin triplet and have 
flavor content {\uu) — \dd))/\/2. Things get more complicated when we consider the mesons rj, 
Tj' and u, (p, however. In the SU{3) flavor-symmetry limit, the pseudoscalar meson is a flavor 
octet and t]' a flavor singlet. However, it is known empirically that S't/(3)-breaking corrections 
to these assignments are large. In the following we shall not rely on SU{3) flavor symmetry, 
but instead introduce another assumption, expected to be accurate at the 10% level. In the 
absence of the axial anomaly, the flavor states \r]q) = {\uu) + \dd))/\/2 and rjs) = |ss) mix 
only through OZI-violating effects, which are known phenomenologically to be small. It is 
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therefore reasonable to assume that the axial anomaly is the only effect that mixes the two 
flavor states j^, . This assumption implies, in particular, that the vector mesons ca and 

0 are pure {\uu) + \dd))/\/2 and |ss) states, respectively, as is indeed the case to very good 
approximation. The anomaly introduces an effective mass term for the system of r] and rj' 
states, which is not diagonal in the flavor basis {\rig), |?7 s)}. Since this is by assumption the 
only mixing effect, one obtains a mixing scheme with a single mixing angle in the flavor basis. 

As explained in js^, the rj and rj' mesons have a leading-twist two-gluon LCDA besides the 
LCDAs corresponding to the quark-anti-quark Fock states r]q and rjs. The two-gluon LCDA 
contributes to the Z —)■ decay amplitudes at order as, through fermion box graphs with 

Zjgg as external particles. A detailed analysis of these decays will be presented elsewhere (33| . 


3 Radiative decays of electroweak gauge bosons 

We now apply our general approach to study the rare, exclusive radiative decays Z —)■ My 
and W —>■ My, where M denotes a pseudoscalar (P) or vector meson {V). The leading-order 
Feynman diagrams contributing to the hrst process were already shown in Figured] We only 
consider cases where the mass of the hnal-state meson is much smaller than the mass of the 
decaying boson. Up to corrections of order {ruM/mz^wY ^^is mass can then be set to zero. 


3.1 Radiative hadronic decays of Z bosons 

We begin our analysis with the decays Z^ —)■ M°y. We hnd that, at leading order in the 
expansion in Aqcd/?7^z, only pseudoscalar or longitudinally polarized vector mesons can be 
produced. The corresponding decay amplitudes can be written in the general form 


iA{Z —)■ My) 


egfhd 

2 cos 9w 




ef 

k ■ q 


pM 


(^z ■ ^ 


* 

7 


q-ezk-e*^\ 
k-q ) 


ipM 
^2 5 


(33) 


where the upper (lower) sign refers to the case where M = P (V||). Here 6\y is the electroweak 
mixing angle. Both the photon and the Z boson are transversely polarized with respect to the 
decay axis. The second term inside the brackets can be written more compactly as and 

below we use this as a short-hand notation. We use a convention where 60123 = — 1- For neutral 
mesons that are eigenstates of the charge-conjugation operation, C invariance implies 

= 0. (34) 


The decay amplitudes are then proportional to the vector product £z x e* of the transversely 
polarized photon and Z boson. However, in new-physics models in which the Z boson has 
flavor-changing neutral-current (FCNC) couplings, Z —)• My decays into mesons that are not 
flavor diagonal (and hence not eigenstates of C) can occur. In this case relation fl3T|) no longer 
holds. In complete generality, the decay rates, summed (averaged) over the polarization states 
of the photon {Z boson), are obtained as 

T{Z M-f) = {\F^\^ + |F^2^r) • (35) 
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Here a = 1/137.036 is the fine-structure constant evaluated at = 0 as appropriate for 
a real photon, and v denotes the Higgs vacuum expectation value, which enters through the 
relation [gj cos6*w)^ = evaluated at /i = mz- This can be solved to give 


V = v{mz) = mz 


sin 6w cos 6 


w 


= 245.36 GeV, 


(36) 


^/Tralmz) 

where we have used a{mz) = 1/127.940 ± 0.014 and 6^ = 0.23126 ± 0.00005, with the 
weak mixing angle determined from the neutral-current couplings of the Z boson evaluated 
at /i = mz ji^. The form factors are given in terms of overlap integrals of calculable 
hard-scattering coefficients with LCDAs. 

Evaluating the diagrams shown in Figures [1] and |2l we hud that the relevant hard-scattering 
coefficients for the decays V —)■ My (with V = Z, W) are given by 


H±{x, mv, P-) = - 

X 


1 h±{x, mv, /i) + 


47r 


where 


h±{x,mv, fi) = (21nx-|-3) ^ 


In 4 

/i 2 


— in] -I- In^ a; — 9 -I- (±1 — 2) 


a; In a; 


(37) 


(38) 


1 — a; 

Our result for h_|_ agrees with a corresponding expression derived in the context of a study of 
meson-photon transition form factors at high performed in 85(]. The expression for is 
new. The relevant convolutions of the hard-scattering coefficients H±{x,mv, with LCDAs 
give rise to the master integrals (we define a^(/i) = 1) 


IM 


fM 

-'it 


(my) = / da;M±(a;,my,/a)0M(a:,/a) = 3 ^(-l)"C'/=^^(my,p)a(f(/a), 


n=0 


(39) 


(my) = / dx H±{1 - x,mv, fi) (j)M{x, ^) = 3C^\mv, fi) (fi), 


n=0 


with 


cy (.nv, ^) = 1 + ey + 0(y), 


An 


\ fi J 


(40) 


The integrals arise from the diagrams shown in Figure [21 in which the photon is attached to 
the quark inside the meson. Diagrams in which the photon is attached to the anti-quark give 
rise to the integrals . In evaluating the integrals we have used the Gegenbauer expansion 
flT3l) . The two types of integrals are related to each other by the fact that the Gegenbauer 
polynomials Cn ' {2x — 1) transform into themselves times a factor (—I)"" under the exchange 
of X O (1 — x). Notice that at tree level the master integrals involve the infinite sums over 
Gegenbauer moments with equal coefficients. Employing a technique explained in Appendix lO 
we have succeeded to derive a closed expression for the one-loop coefficients c^\mv/l^)- It 
reads 


44,()^). 


_(n + l)(n 2) 

4(M„+i 


-|- 3 


In 4 


in 


AH, 


1)±1 


( 41 ) 


n +1 


+ 


(n-I-l)(n-I-2) {n + iy{n + 2y 


9. 


17 




















0.60 r 


3. 0.58 


"5 


0.56 


0.54- 


O 


0.52 


0.50 ^ 


50 


75 100 125 150 175 

ji [GeV] 



H [GeV] 


Figure 4: Scale dependence of the combinations C^\mzi jj) {jj)/{fio) for the 
hrst two Gegenbauer moments (n = 1,2). The red dashed lines show the results at 
leading-order, where C^\mz, = 1- The blue and yellow lines show the results at 
next-to-leading order obtained when the one-loop expressions in (HUjl are used. 


Note that = —5 and Cq~^ = —4 are pure numbers. Using the evolution equations (12^ and 
the explicit expressions for the one-loop anomalous dimensions given in fl30p . it is straight¬ 
forward to check that the master integrals in fl3^ are independent of the factorization scale 
p. Indeed, the coefficient of the logarithm in flTTl) is equal to —7„/(2Gi?). Note also that the 
imaginary parts associated with the logarithm do not contribute to the decay rates at 0(as). 

In Figure IH we study the scale dependence of individual terms in the sums over Gegen¬ 
bauer moments in fl5^ at leading (dashed red lines) and next-to-leading order (solid lines) in 
perturbation theory. At leading order the /i dependence of the Gegenbauer moments, shown 
explicitly in is left uncompensated, and hence a signihcant scale dependence arises. At 
next-to-leading order this dependence is compensated by the logarithmic terms contained in 
the one-loop corrections (ITT]) to the hard-scattering coefficients Ci'^\mz, ju) (blue lines) and 
Ci~^(mz, ju) (orange lines). The resulting next-to-leading order curves exhibit excellent sta¬ 
bility under variations of the factorization scale in the interval m^/2 < /i < 2mz- 

In terms of the master integrals dehned in fl39l) . the form factors are given by 


O °° 

^ [I^(mz) + If(mz)] = Qm'^ C^^^mz,/^) a^(p), 

n=0 

O' °° 

[l-{mz) - I^{mz)\ = -Q'm a^+i(h), 

n=0 


(42) 


where 

Qp = '^ 6CgQgVg, Qv = ^QCgQqaq, (43) 

<? <? 

and the coefficients and related to Qm by exchanging Vg ^ Qg. Gorrections to the results 
flT2]) arise only at twist-4 level and are suppressed by (Agco/^TT-z)^ or {ruM■ They are 
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Meson P 

Qp = 'EqQc^Qq'Vq 

Meson V 

QV = Y.q^C^Qq(^q 

7r° 

^(l-4sin2 9w) 


1 

272 

Vq 


UJ 

3 

2^/2 

Vs, Vb 

1 — 1 sin^ 6w 

0, T 

1 

2 

Vc 

1 — 1 sin^ 9w 

J/0 

1 


Table 3: Coefficients Qm for the ground-state neutral pseudoscalar and vector mesons. 


phenomenologically irrelevant. In the above expressions Qq denotes the electric charge of a 
quarks in units of e, while Vq = — sin^ 6w Qq and Qq = \T^ (not to be confused with 

the Gegenbauer moments) are its vector and axial-vector couplings to the Z boson. Our 
hnding that the form factors for pseudoscalar and vector mesons in fl42p have exactly the 
same structure crucially relies on a mathematically consistent treatment of 75 , see Section 1^751 
At tree level = 1, and hence the form factor {F^) is proportional to the inhnite sum 
of all even (odd) Gegenbauer moments of the meson M. Gharge-conjugation invariance implies 
that the LGDAs of a flavor-diagonal neutral mesons are symmetric under the exchange of x and 
(1 — x), and hence for these mesons the odd Gegenbauer moments vanish. This leads to 
relation fl5T)) . The non-zero form factor involves the inhnite sum over the even Gegenbauer 
moments times some havor-dependent coefficients Qm, which we collect in Table |3l 

Explicit predictions for the leading-twist LGDAs derived by means of non-perturbative 
methods are typically obtained at a low hadronic scale /xo ~ 1 GeV. When these predictions are 
used in fH 2 |) . the expressions for the radiative corrections involve large logarithms ln(m|//iQ) 

9, which must be resummed to all orders in perturbation theory to obtain reliable predictions. 
This resummation is most readily performed by evaluating the result 042 p at the scale /i = mz 
(or any other scale of the same order), in which case we obtain 


ReFj^ = Qm 0.94 1.05 a“(mz) 1.15 a^(m^) -h 1.22 a“(m^) ... 




Ml 


.Ml 


= Q 


M 


0.94 + 0.41 af (/io) + 0.29 af (/xq) + 0.23 af (/xq) + 


(44) 


We use the three-loop expression for the running coupling as provided by the RunDec pro¬ 
gram j^, normalized to as{mz) = 0.1185 ± 0.0006 jdsj and with heavy-quark thresholds at 
mb{mb) = 4.163 GeV and mc{rnc) = 1.279 GeV j^. The Gegenbauer moments at the high 
scale /X = mz in the hrst line can be related to hadronic input parameters calculated at the 
low scale /xq = 1 GeV using the relations fl3Tp . In this process the coefficients of the higher 
moments get successively smaller. 

Decays into a transversely polarized vector meson are only allowed at twist-3 order. This 
presents us with an opportunity to study the structure of power corrections with a specific 
test case. We adopt the approximation where three-particle LGDAs are neglected. We then 
evaluate the diagrams in Figured] using the projector for a transversely polarized vector meson 
given in Appendix]^ The decay amplitude can be decomposed in a form analogous to fl5^ . 
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such that 


iA{Z\\ -)■ Kl 7 ) 


egfv my 
2 cos 9w fnz 




k^q 


1 /^* 0:-*/3 

k ■ q 


F,^-s 
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(45) 


The Z boson must be longitudinally polarized, and its polarization vector can be written 
as = (g — kY/mz- The extra factor of my jraz compared with fl33|l makes the power 
suppression of these amplitudes explicit. The corresponding decay rate, summed (averaged) 
over the polarizations of the final-state (initial-state) particles, is given by 


r(z ^ 147) 


amzfy ml /, ^,2 ^ |^^| 2 \ 
m| V' ^ ' I 2 I y 


(46) 


The general expressions for the form factors FY in terms of overlap integrals over the various 
twist-3 LCDAs appearing in the projector flA.4p for a transversely polarized vector meson are 
given in Appendix They can be simplified a lot by using relations implied by the equations 
of motion in the limit where three-particle LCDAs are neglected. Assuming for simplicity that 
quark-mass effects can be neglected, we obtain 


= 


Q 


V 


dx 


In a; ln(l — x) 


1 — X 


+ 


X 


(j)y{x,fi) + 3-particle LCDAs . 


(47) 


An analogous expression with Qy replaced by Q'y and a relative minus sign between the 
two terms inside the parenthesis holds for F^- Since the leading-twist LCDAs of neutral 
mesons are symmetric in a; 77 (1 — x), it follows that F^- = 0, and hence once again only the 
term involving the Levi-Civita tensor in fl4^ contributes to the decay amplitude. Using the 
Gegenbauer expansion of the LCDA 0y(x,/a) in flT^ . we obtain 


fY = q 


V 




4 

^2n 




^ (n + l)( 2 n + l) 


(48) 


Since we have not evaluated radiative corrections to the form factors, we do not control the 
scale dependence of the Gegenbauer moments a)(. However, it is clear that in order to avoid 
large logarithms we should again set /r ~ in the final result. Note that the result for F^ 
is similar to that for F^ in 042 p , but the coefficients of higher Gegenbauer moments are more 
strongly suppressed. For a rough estimate, we may assume that the Gegenbauer moments 
only have a minor impact on the final results (i.e. al{mz) S> 1 ), in which case it follows that 
the rates for Z —)■ U 7 decays with transversely and longitudinally polarized vector mesons are 
related by 


T(Z ^ V7) ,, ml 
T(Z ^ V)| 7 ) ~ ml ■ 


(49) 


This ratio is of order 10 The outcome of this discussion is that power corrections in the 
expansion in Aqcy,/mz are completely negligible for phenomenological applications. 
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Figure 5: Non-local (left and center) and local (right) contributions to the W~^ —)■ 
decay amplitudes. 


3.2 Radiative hadronic decays of W bosons 


The exclusive radiative decays W~^ —)■ M+y are, at hrst sight, very similar to the decays 
—)■ M^'j. Indeed, the contributions from the hrst two diagrams shown in Figure [5] can be 
obtained from the corresponding contributions to the Z-boson decay amplitudes by means of 
simple substitutions. The charged currents are havor non-diagonal, and hence the hnal-state 
meson M~^ has a dehnite havor structure described by a wave function \uidj). Note that now 
diherent electric-charge factors arise, depending on whether the photon is attached to the up- 
type quark or down-type anti-quark. The charged currents are purely left-handed, and hence 
we must replace 


Vq, ttq —?• 


COS 9w 
2y/2 




(50) 


in the equations of the previous section. However, a careful analysis shows that the hrst two 
diagrams in Figure E] give rise to an extra contribution with a diherent tensor structure. It 
reads 


iAA{W+ M+y) 


egfiM 


^ij i^Qu 


Qd) 


k-e*^q-ew 
k ■ q 


(51) 


where the upper (lower) sign refers to the case of a pseudoscalar (longitudinally polarized 
vector) meson in the hnal state. Note that this contribution is independent of the LCDA of 
the hnal-state meson. It vanishes for an on-shell (transverse) photon, but is not compatible 
with f/(l)em gauge invariance. 

Since the W boson has a direct coupling to the photon, an extra contribution to the 
—)■ M+y decay amplitudes exists, which arises from the third diagram in Figure 13 in 
which the hnal-state meson is produced by the conversion of an oh-shell W boson. This graph 
has no analog in the Z-boson case. The corresponding contribution to the decay amplitude 
involves the meson matrix element of a local current, which to all orders in QCD is given in 
terms of a meson decay constant. We hnd 


^ P+y) 


^Zllocal(b"+ ^ V^+7) 


egfp 

2 V 2 


Vij Ew 



( 52 ) 
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where we keep the exact dependence on the vector-meson mass tum for the time being. The 
second relation can be simplihed by considering the cases of longitndinal and transverse po¬ 
larization separately. The polarization vector for a longitndinally polarized vector meson can 
be decomposed as 


£ 


IIm 

V 


1 

my 



2my 

^w-^v 



(53) 


which satishes the conditions k ■ = Q and {sy)‘^ = —1. The polarization vector for a 

transversely polarized vector meson is dehned snch that k ■ Sy = q ■ Sy = We then obtain 


Ulocal(hh+ ^ K+7) 


*^local(fh|^ ^ fo/7) 


2^/2 


By/ -8*^ + 0 



<^9fv y mv _L* 

2 V 2 ^ 



(54) 


The second amplitnde is non-zero only if the W boson is longitndinally polarized, and we have 

in the hnal resnlt. The 


my 


used a decomposition analogous to fl5^ to replace —2q ■ bw = 
local amplitudes for M = P,V» are such that they combine with the extra term in fl5T|) to give 


a gauge-invariant result proportional to b^ ■ b-^* [24, [25 

It follows from this discussion that, in analogy with fl5^ . the leading-power amplitudes for 
the decays —)■ M+y can be written in the general form 


a(iv+ ^ m+7) = r., ■ d’ 7') ■ (ss) 

Summing (averaging) over the polarization states of the photon (W boson), we obtain the 
corresponding decay rates 

r(W+ ^ m+ 7) = Ifo.f (ifof I' + • (56) 

In close analogy with fl42p . we hnd that the form factors are given by 


= QuI+ijnw) + Qdl+ijnw) = ^ a^(/i) - , 9) a^+i(h) , 

Ff = -2 (Q, - Qd) + Qu - Qd I^{mw) (57) 

CO 

= -2 + ^ ^3Ciz\rnw,^i) a^+i(h) • 

n=0 

The contribution —2 to arises from the local contribution in fl52l) . Corrections to these 
results are suppressed by (Aqcd/'"2-w)^ or The corresponding amplitudes for the 

decays W~ —)■ M “7 are obtained by replacing Idj —t Id* in fl5^ and by replacing the charge 
factors Qu 77 Qd in fl57)) . In addition, one must take into account that the odd Gegenbauer 
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moments of the meson M~ have the opposite sign as those of . This can be accounted for 
by replacing . As a result, the form factor remains invariant, while the form 

factor changes sign. The decay rate in fl56ll stays invariant under these replacements. At 
low values of the factoriszation scale p, the Wilson coefficients in ll57|l contain 

large logarithms of the form ln(m^//i^), which can be resummed to all orders in perturbation 
theory by evaluating the scale-invariant quantities in flFr|) at the scale fi = mw We obtain 

Re F^ = 0.94 — 2.98 + 1-05 a^{miv) — 3.31 a^{miv) + 1.15 af{mw) =F • • • 

= 0.94 - 1.65 af (/io) + 0.42 af (/xq) - 1.03 af (/xq) + 0.30 af (/io) T • • • , 

(58) 

Re F^ = 0.85 — 1.00 {mw) + 3.16 af {mw) — l-H af {mw) + 3.45 af {mw) T • • • 

= 0.85 - 0.55 af (/io) + 1.25 af (/xq) - 0.34 af (/xq) + 0.89 af (/xq) T • • • • 

Decays into a transversely polarized vector meson are once again only allowed at twist-3 
order. We decompose the decay amplitude in a form analogous to fl45D . such that 


iA{Wf 1/+7) 


eg'/v my 
Af2 mw 






k ■ q 


Ff - e 


V 



(59) 


The corresponding decay rate, summed (averaged) over the polarizations of the hnal-state 
(initial-state) particles, is given by 
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m 
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7 + 7 


(60) 


The form factors can be calculated in analogy with the discussion in the previous section. The 
hnal results are 


= —2 / dx 


Qu Inx Qdln(l - x) 
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0 v()r, /x), 


Ff = -A{Q^-Qd)-2 / dx 


Qu^^x Qd\a{l-x) 


1 — X 
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(61) 


0y(a;,/x), 


where for simplicity we neglect contributions proportional to the quark masses. The local 
contribution in (154)) adds —2{Qu — Qd) to Ff. When expressed in terms of Gegenbauer 
moments, these results take the form 




®2n(h') 


-E 


3r7^ 

'J“ 2 n+l 


(h) 


^(^+l)(2n + l) ^(n + l)(27x + 3)’ 


^2^ = -1 + E 


3®2n(h') 


(n + l)(2n + l) ^ (n + l)(2n + 3) ■ 

= 1 ^ ^ ^ ' n =0 ^ ^ ^ ^ 
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a. 


V 

2n+l 


(h) 


In the limit where the Gegenbauer moments are neglected, we hnd in analogy with 
the ratio 


r(W+ ^ R+7) ml 

T{W+ f+7) ~ ^ 


(62) 

that 

( 63 ) 


is strongly suppressed. 


23 

















'X/X/V/Xz 7 


7 



Figure 6: One-loop triangle graphs giving rise to the axial anomaly. The dot represents 
the axial-vector current Aj^. 


3.3 Absence of enhanced contributions from the axial anomaly 

In has been suggested in the literature that the decay amplitudes for Z —)■ Py and W 
Py receive a very large enhancement due to an analog of the axial anomaly, which gives a 
contribution to the 7r° —)■ yy decay amplitude that does not vanish in the chiral 

limit [3, [ 2 ^. We now explain why such a contribution does not exist in our case. 

We consider the axial current = q'j^'y^q and evaluate the triangle diagrams shown in 
Figure [HI where instead of two photons we take the external particles to be a photon and a Z 
boson. When taking the divergence of the current and considering the chiral limit —)■ 0, 
we hnd that the loop diagrams vanish if one naively assumes that ys anti-commutes with 
all Dirac matrices y^. However, a more careful regularization prescription adopting the HV 
scheme shows, like in the case of two external photons, that a hnite remainder exists. It 
corresponds to a local operator for two gauge bosons. In operator language, we hnd that 

a,A!; = 2tm, m ^ W (qI f ""Z"'’ + (64) 

^ 47r \^ smOwcosOw J 


where for completeness we have included the terms proportional to the quark mass on the 
right-hand side. We do not show contributions involving two electroweak gauge bosons {ZZ 
or W^W~) or two gluons, since they are irrelevant to our discussion. After electroweak 
symmetry breaking, the heavy Z boson acts like an external source, which is invariant under 
P(l)em gauge transformations. We see that anomalous contributions to the divergence of 
the axial-vector current not only involve the photon held, but also the Z boson. Indeed, 
such anomalous terms even arise for the charged currents Ah = d^y^ysUi. Evaluating the 
corresponding triangle graphs with an external ylF state, we obtain 

= i{mui + md^)dj'y 5 Ui + ie dj4lbUi - e^^ai3 P^‘'IF+“^ -F ... , (65) 

47r 3v2sm6'vK 


where we omit a contribution involving two electroweak gauge bosons (WZ). 

Does the existence of the anomalous di-boson terms in the above relations imply that there 
exist enhanced contributions to the Z —)■ Py and W —)■ Py decay amplitudes, in analogy to the 


famous case of the vr^ —)■ yy amplitude? This possibility was suggested in |28|, [29 


where the 

authors speculated about a huge enhancement of the rates for the radiative decays Z —)■ vr^y 
and W~^ —)■ Dg'j (see also the recent paper [l^, which considers the decay W~^ 7r“y). We 
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decay of a Z boson into a photon and an axial current Right: A hypothetical 
anomaly-mediated contribution to the —)■ Pj decay amplitude. The crossed circle 
represents the decay amplitude in flBBl) : the double line shows the meson propagator. 


now demonstrate that such an enhanced contribution does not exist, focussing for concreteness 
on the case of a neutral pseudoscalar meson meson P. Let us parameterize a hypothetical 
anomalous contribution to the Z —)■ Pj decay amplitude in the form 

iA(Z Pj) = iA , (66) 

where the general structure of the amplitude is consistent with fl33|l . Let us furthermore 
parameterize the amplitude coupling an initial-state Z boson to the axial current A^ and a 
photon as 

= (67) 

At lowest order, this amplitude is obtained from the diagram shown on the left in Figure [71 
Inserting a complete set of hadron states that can be interpolated by the axial current A^, and 
summing over quark flavors, we find that the amplitude iAi^'^^{k,q) contains the following 
contribution from the the single-hadron state P: 

q)^'^Cg (-i/p/c^) ^ k^ep^afi k^q^ . (68) 

q P q 


This can be read off from the graph shown on the right in the figure. In the last step we have 
taken the chiral limit, in which the meson becomes massless. The key feature of the anomalous 
contribution would be that it exhibits a 1/fc^ pole in this limit. We have calculated the one- 
loop contributions to the amplitude by evaluating the loop diagrams in Figured We 

hud that the contribution associated with the tensor structure shown in fl68|) is proportional 
to the expression (with p = k + q) 


F(l + e)F^(l-e) 

F(2 — 2e) (p2 _ 


, - (-kT‘ . (-kp--- (-pp-‘ 

■t 1 

e 1 — e 


(69) 


where + iO and k"^ = k"^ + iO. For the case of the 7r° —)■ 77 amplitude = 0 for the 

external photon, and in the limit e —)■ 0 one obtains —1/fc^, which indeed exhibits a pole. Form 
the residue of this pole one can derive the anomaly-mediated tt^ —)■ 77 decay amplitude. For 
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Decay mode 


Branching ratio 


asymptotic LO 


—>■ 7r°7 

ZO ^ pOy 

07 

ZO ^ J/07 
ZO ^ T(1S)7 
ZO ^ T(4S)7 


(9.80 ± 0.03/ ± 0.61,^ ± 0.82„J • lO’^^ 

(4.19± 0.16/ ± 0.24„, ± 0.37,J • lO’^ 
(2.821°;°!^ ± 0.15/ ± 0.28„, ± 0.25,J ■ lO’^ 
(1.04 + °;°! ^ ± 0.05/ ± 0.07„, ± 0.09,J ■ lO’^ 


-0.06 


0.02 u.uo/ _i_ u.u I a2 


7.71 14.67 

3.63 5.68 

2.48 3.76 

0.86 1.49 

10.48 6.55 

7.55 4.11 

1.71 0.93 

13.96 7.59 


ZO ^ T(nS )7 



Table 4: Predicted branching fractions for varions Z —)■ M'j decays, inclnding error 
estimates dne to scale dependence (snbscript “/i”) and the nncertainties in the meson 
decay constants (“/”), the Gegenbaner moments of light mesons (“a„”), and the width 
parameters of heavy mesons (“cr”). See text for fnrther explanations. 


onr case, on the other hand, = m\ is eqnal to the mass of the decaying heavy gange boson, 
in which case the above expression does not exhibit a 1/fc^ pole, bnt is instead proportional 
to l/ml. Hence we conclnde that H = 0 in fl65]) . Note that in the limit —)■ 0 one obtains 


from fl69D 



(70) 


which is precisely of the form of onr (bare) hard-scattering coefficients. 

3.4 Phenomenological results 

We are now ready to present detailed nnmerical predictions for the various radiative decay 
modes. We start with the decays of the Z boson, using relation (I5B]) . Besides the input 
parameters already mentioned, we need the Z-boson mass mz = (91.1876 ± 0.0021) GeV and 
total width Pz = (2.4955 ±0.0009) GeV ji^. When squaring the decay amplitudes, we expand 
the resulting expressions consistently to hrst order in ag. The imaginary parts of the form 
factors in flT2l) do not enter at this order. Our results are presented in Table HI Signihcant 
uncertainties in our predictions arise from the hadronic input parameters, in particular the 
meson decay constants (see Appendix [B]) and the various Gegenbaner moments. Their impact 
is explicitly shown in the table. Our error budget also includes a perturbative uncertainty, 
which we estimate by varying the factorization scale by a factor of 2 about the default value 
p = mz- All other uncertainties, such as those in the values of Standard Model parameters, 
are negligible. Note also that power corrections from higher-twist LGDAs are bound to be 
negligibly small, since they scale like (Aqcd/'^z)^ for light mesons and at most like (mM/mz)"^ 
for heavy ones. The predicted branching fractions range from about 10“^^ for Z° —)■ to 
about 10“^ for Z° —)■ J/'ip'j. In the last row, the symbol T{nS) means that we sum over 
the hrst three T states (n = 1,2,3). Strong, mode-specihc differences arise foremost from the 
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Figure 8 : Examples of QED (left) and electroweak radiative corrections (center and 
right) to the Z —)■ M 7 decay amplitudes. The last two diagrams can give rise to 
flavor-violating decays in the Standard Model. 


relevant flavor-dependent coefficients in Table [3l as well as from differences in the values of the 
decay constants. The combined uncertainties in the predictions for the branching fractions are 
typically of order 10% and are dominated by the uncertainties in the shapes of the LCDAs. 
The only exception are the decays —)■ T 7 , for which the relevant hadronic overlap integral 
is constrained by the model-independent relation ffTOjl . 

In our analysis we neglect two-loop QCD corrections, whose effects should be covered by the 
error we estimate from scale variations, and one-loop QED or electroweak radiative corrections, 
a few examples of which are shown in Figure |H1 Their impact should be much smaller than 
the theoretical uncertainties inherent in our predictions. Consider, as a concrete example, the 
contribution of the hrst diagram, which only contributes to the Z —)■ I /7 amplitudes. Notice 
that the photon propagator l//c^ with is cancelled, because the Z —)■ 77 * amplitude 


vanishes if both photons are on-shell 88 |. As a result, there is no enhancement factor and the 


diagram is suppressed, compared with the leading contributions shown in Figure [H by a factor 
a/vr ~ 2 ■ 10“^. This naive estimate is conhrmed by the result of a detailed calculation of this 
contribution to the Z J/4’1 decay amplitude, which found that its effect leads to a reduction 
of the leading contribution by 0.2%, corresponding to a 0.4% correction of the branching ratio 
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In the same paper, the authors have presented predictions for three of the Z —)■ I /7 decay 
modes along with theoretical error estimates. They are Br(Z (p'j) = (11.7 ± 0.8) ■ 10“®, 
Bt{Z J/'ip-i) = (9.96 ± 1.86) • 10-®, and Br(Z ^ T(15)7) = (4.93 ± 0.51) ■ 10"^ The 
last two branching ratios are consistent with our Endings within errors. Note that in the 
NRQCD approach adopted by these authors the decay constants of the heavy quarkonia 
are themselves derived from an expansion about the non-relativistic limit. This introduces 
additional uncertainties, which can be avoided if the decay constants are extracted from data, 
as discussed in Appendix [Bl The analysis of the decay Z^ —)■ 07 presented in 
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uses an 

approach similar to ours but only includes the leading logarithmic evolution eEects from the 
hadronic scale po = 1 GeV to the high scale /i ~ mz- Their result is consistent with ours 
but has a smaller uncertainty. The non-logarithmic 0{as) corrections included here for the 
first time reduce the branching ratio by a significant amount. We also find that the present 
ignorance about the precise shape of the 0-meson LCDA gives rise to a larger uncertainty. 

We now proceed to present our predictions for exclusive radiative decays of W bosons. 
In this case we need the input parameters mw = (80.385 ± 0.015) GeV and Tw = (2.0897 ± 
0.0008) GeV, as well as the relevant entries of the quark mixing matrix, which are \Vud\ = 
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Decay mode 

Branching ratio 

asymptotic 

LO 


—>■ 7r=*=y 

(4.00 

± 0.01/ ± 0.49^2 ± 

0.66,J . 

10- 

-9 

2.45 

8.09 


^ P^y 

(8.74 

± 0.33/ ± 1.02^2 ± 

1.57„J . 

10- 

-9 

6.48 

15.12 


^ A±y 

(3.25):»S,±0. 

33/ ± 0.24„, ± 0.38 

± 0.51 

114 ) 

10-10 

1.88 

6.38 


A*±y 

(4.781^4. ±0. 

28/ ±0.39ai ±0.66 

a2 ± 0.80 

114 ) 

10-10 

3.18 

8.47 


-t Ds'y 

(3.66 

± 0.12ckm ± 0.13/ 

+ 1.47 1 
-0.82 rrj ' 

10- 

-8 

0.98 

8.59 


—)■ D^y 

l^i.OO _o.02 M 

± O.lOcKM ± 0.07/ 

+ 0.50 I 
-0.30 o-J ■ 

10- 

-9 

0.32 

3.42 


^ R±y 

71 ccr+O.OO 
-0.03 /I 

±0.37ckm±0.15/ 

+ 0.68 \ 
-0.45 o-j • 

10- 

12 

0.09 

6.44 


Table 5: Predicted branching fractions for varions W —>■ M 7 decays, inclnding error 
estimates dne to scale dependence and the nncertainties in the CKM matrix elements, 
the meson decay constants and the LCDAs. The notation is the same as in Table 01 
See text for farther explanations. 


0.97425 ±0.00022, |14.| = 0.2253 ± 0.0008, |14s| = 0.986 ±0.016, \Vcd\ = 0.225 ±0.008, iKbl = 
(41.1 ±1.3)-10“^, and |14b| = (4.13±0.49)-10“^ ji^. Starting from relation fl56|l . we obtain the 
resnlts shown in Table O In this case the pattern of the different decay modes reflects mainly 
the pattern of the relevant CKM matrix elements, and to a lesser extent the differences in the 
decay constants. The Cabibbo-allowed decays W —)■ Try, py, and H^y have branching fractions 
of order few times 10 “® to few times 10 “®, where decays into heavy mesons are enhanced 
dne to the strnctnre of the relevant overlap integral in ffT 8 |) . The Cabibbo-snppressed modes 
W y and the strongly CKM-snppressed decay W Bj have correspondingly smaller 

branching ratios. The nncertainties inherited from CKM elements are shown where they are 
signihcant. In a recent paper, the —>■ vr^y branching ratio was estimated to be 0.64 • 10“® 
1^, which is abont 6.3 times smaller than the valne we obtain (see below). 

In the last two colnmns in Tables 0] and [5] we show different approximations to onr resnlts. 
The hrst one (labelled “asymptotic”) gives the central valnes of the branching ratios (in the 
appropriate nnits) obtained if the asymptotic form 6 a:(l — x) of the meson LCDA is employed. 
As we have explained, RG evolntion effects from the low hadronic scale po = 1 GeV np to 
the electroweak scale have the effect of strongly snppressing the contribntions from higher 
Gegenbaner moments. Indeed, we observe that nsing the asymptotic form provides reason¬ 
able approximations in most cases (especially for the Z —)■ My modes). The corresponding 
expressions for the decay rates read 


r(Z° ^ M°y)| 
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^ arnzfl 2 
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where Vij is the relevant GKM matrix element for the prodnction of the charged meson M~^. 
The dominant corrections to the Z —)■ My branching fractions arise from the second Gegen¬ 
baner moment a^, which is positive for light mesons and negative for heavy qnarkonia. The 

















Decay mode 

Branching ratio 

SM background 

Z^ ipOy 

^ D°y 
^ H°y 

[(7.70 ± 0.83) \v,X + (0.01 ± 0.01) |a.i|2] . 10-« 
[(5.30 ISS) |t>o.P + (0.62y:g) |o„p] . 10-' 
[(2.08^112) |t.Mp + (0.77iS) |o„p] ■ 10-' 
[(2.64111:1) + (0.87lll:g) |o„p] . 10-' 

-4^ “ ~ 2 • 10-3 

Sin U\Y 77 

-4^ “ ~ 2 ■ 10-3 

sin 77 

~ 8 • 10-3 
-4^ “ ~ 4 ■ 10 -^ 

sm dy/ 77 


Table 6 : Branching fractions for FCNC transitions Z —)■ M 7 , which conld arise from 
physics beyond the Standard Model. The different theoretical nncertainties have been 
added in quadratnre. The last colnmn shows onr estimates for the irredncible Standard 
Model backgronnd np to which one can probe the flavor-changing conplings and 
aij. Here A ~ 0.2 is the Wolfenstein parameter. 


dominant corrections to the W —)■ Mj branching fractions with kaons, D mesons or B mesons 
in the hnal state arise from the hrst Gegenbaner moment . It gives a large positive contri¬ 
bution in all cases. In the case of heavy mesons this effect is particularly pronounced. The 
approximate results in (ITT]) are fully consistent with corresponding (tree-level) expressions 


derived in [2^. The result for the —)■ decay rate derived in [2^ is lower than ours by a 

factor 4/9, and the formula for the —?• vr^y decay rate derived in jl^ differs from 0711) by 

a factor 2/9. The origin of the discrepancy is related to the fact that the theoretical approach 
used in these papers is based on an expansion in a parameter uq = 2, and the numerical 
estimates are obtained by keeping only the leading term in the expansion - a fact that was 
admitted in these papers. In Appendix [D] we trace the source of the discrepancy in more 
detail. In the last column in Tables H] and 0 (labelled “LO”) we present the branching ratios 
one would obtain at tree level using the model predictions for the LCDAs at the low scale /iq. 
In this approximation the one-loop QCD corrections, which contain large logarithms of the 
form tts ln(m|^/pg), are omitted. For most decays the corresponding results overshoot the 
values obtained at next-to-leading order by signihcant amounts; only for decays into heavy 
quarkonia they underestimate the branching fractions. 

Future precision measurements of the exclusive radiative decays Z —)■ My would serve 
as powerful tests of the Standard Model and of the framework of QCD factorization. The 
branching ratios for decays into vector mesons shown in Table [Hare proportional to |agp, where 
ttg denote the axial-vector couplings of the Z boson to the various quarks. The couplings \ab\ 


and lad have been measured at LEP with 1% accuracy [89|], but no similarly accurate direct 
measurements of the couplings to light quarks are available. Given the theoretical precision 
of our predictions, it would be possible to determine these couplings with about 6 % accuracy. 
The decays Z —)■ My can in principle also be used to search for non-standard FCNC couplings 
of the Z boson. If such couplings exist, then the diagrams shown in Figured] can lead to hnal- 
state mesons of mixed flavor, such as D^, and Bg. It is straightforward to calculate 
the corresponding decay rates in our approach, starting from the general relations fl35p and 
m- We parameterize the non-standard vector and axial-vector couplings of the Z boson 
by Vij and aij, respectively, where i,j are the quark flavors of the hnal-state meson. Our 
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|Re[(usd ± asd)"^] 1 

< 2.9 ■ 10-® 

lRe[(usd)2 - (usd)^] 1 

< 3.0 ■ 10-1° 

|lm[(u^d ± ttsd^] 1 

< 1.0 ■ 10-1° 

|lm[(u^d)2 - (usd)^] 1 

< 4.3 ■ 10-13 

1 {^^CU i ^cu) 1 

< 2.2 • 10-® 

\{VcuY - (Ocn)^! 

< 1.5-10-8 

\{vm ± Om)^! 

< 4.3 • 10-® 

livbd)"^ - {abdf\ 

< 8.2 • 10-° 

\{Vbs ± flbs)^ 

< 5.5-10-^ 

(^^6s)^ - (Ofe)^ 

< 1.4-10-^ 


Table 7: Indirect constraints on the flavor-changing Z-boson couplings 
95% confldence level) derived from neutral-meson mixing 90-^. 


and ctj^j (at 


predictions for the corresponding branching fractions are given in Table El At higher order 
some Standard Model background to these searches exists, since electroweak loop graphs such 
as those shown in the last two diagrams in Figure | 8 ] can give rise to flavor-changing transitions. 
Naive dimensional analysis shows that the contributions of these diagrams, relative to the 
contributions from the graphs in Figure El (in units of the new-physics couplings Vij and aij), 
scale like (a/Tr) \ VikV^j\/ siv? 6w-i where k can be any one of the three possible generation 
indices. The relevant loop functions depend on the dimensionless ratios m^/m| and 
which are either of 0{1) or can be set to zero. Consequently one can only probe the new- 
physics couplings Vij and up to some irreducible Standard Model background, which is 
estimated in the last column in Table El 

Possible FCNC couplings of the Z boson are heavily constrained by precision flavor physics, 
in particular by bounds on the AF = 2 mixing amplitudes. It is a straightforward exercise to 
match our parameters Vij and ajj onto the Wilson coefficients in the general effective AF = 2 
Hamiltonian as defined, e.g., in |90|. We obtain 


Ci = 


AGf 


I 

+ (^ij) 


Ci = 


AGf 


V 2 ^ ^ ^/2 

All other coefficients are zero at tree level, 
updated results reported in 


Vi4 Cl' 


VJ 


G,= 


AGf 


(4 - al) . (72) 


Using the bounds compiled in ^ as well as 
9ll. l92| . we find the upper bounds on various combinations of Vij 


and ttij parameters shown in Table [71 The strongest bounds exist for the coefficients C 5 of 
mixed-chirality operators (right column). They can be avoided by assuming that Vij = ±ajj, 
such that the flavor-changing couplings are either purely left-handed or purely right-handed. 
Under this assumption one finds from the table that < 8.5 ■ 10“®, \vcu\ < 7.4 • 10“®, 
|um| < 1.0 • 10“"^ and |u 6 s| < 3.7 • 10“^, and the same bounds apply to \aij\. If these indirect 
bounds are used, then the branching fraction shown in Table El are predicted to be at most 
a few times 10 “^^ (a few times for the case of —)■ Bg'y), meaning that they will be 

unobservable at the LHC and all currently discussed future facilities. We And it nevertheless 
worthwhile to illustrate the general idea of such new-physics searches. First of all, it should be 
emphasized that the indirect bounds derived from K—K, D—D and Bd,s~Bd,s mixing are to 
some extent model dependent, since one cannot tell whether the flavor violation originates from 
the couplings of the Z boson or from some other new particle. It is conceivable that in some 
(admittedly fine-tuned) models flavor-violating couplings of the Z boson can be compensated 
by the effects of some other, heavy boson. Also, in deriving the bounds on a particular Wilson 
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coefficient Q one assumes that a single new-physics operator is present at a time and sets the 
coefficients of all other operators to zero. The method presented here, on the other hand, is 
unique in that it allows one (in principle) to probe for flavor-changing couplings of the Z boson 
directly and in a model-independent way, based on tree-level couplings of an on-shell particle. 
It should thus be seen as a complementary way to search for such effects. This method can 
also be generalized to the interesting case of flavor-changing exclusive Higgs-boson decays 


for which the corresponding indirect bounds have been studied in [93 


4 Weak radiative hadronic decays M^W 


Exclusive decays of a Z boson into a W boson and a single meson M are kinematically allowed 
as long as the hnal-state meson is lighter than the mass difference mz — mw — 10.8 GeV. 
While similar at hrst sight to the radiative Z-boson decays studied in Section [21 these decays 
are nevertheless interesting for several reasons. Unlike the photon, the hnal-state W boson 
can be longitudinally polarized, and hence several different helicity amplitudes contribute 
to the decay. Also, the trilinear ZWW coupling in the Standard Model gives rise to an 
additional contribution to the decay amplitude, in which the hnal-state meson is produced 
via the conversion of a lU boson. This term is analogous to the “local” contribution we 
encountered in our study of the radiative decays of W bosons. Indeed, the leading-order 
Feynman graphs contributing to the Z° —>■ M^W~ decay amplitudes, shown in Figure [9l are 
analogous to those in Figure [5l Finally, and most interestingly, the decays Z —)■ MW oher 
an opportunity to test the QCD factorization approach at a scale signihcantly lower than the 
Z-boson mass. A factorization theorem of the form (|2]) can only be derived if the momentum 
of the hnal-state meson in the rest frame of the decaying particle is much larger than its mass, 
since only then the constituents of this meson can be described in terms of collinear quark 
and gluon helds in SCET. The relevant condition is 

AAliAllA) » , where \(x,v,z) = - y - zf - iyz. (73) 

2mz 

2 2 

This condition is satished as long as tum -C ~ 10.2 GeV. We can thus use the fac¬ 

torization approach to calculate the branching fractions for Z —)• MW decays with a light, 
strange or charm meson in the hnal state, but not with a B meson. It also follows that the 
non-perturbative corrections to the factorization formula are organized in powers of () 
rather than (Aqcd/'^z)^- 

The local contribution to the decay amplitudes involves the meson matrix elements of local 
currents. In close analogy with ([52]), we hnd 
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Figure 9: Non-local (left and center) and local (right) contributions to the Z M^W 
decay amplitudes. 


These results are exact even as far as the dependence on the vector-meson mass is concerned. 
The second relation can be simplified by considering the cases of longitudinal and transverse 
polarization of the vector meson separately. The longitudinal polarization vector can be 
written as 

jl^_ 1 ml - mlr - ml ( 

~ \/2 2 2 \ \ ^ 
mv A(m|,m^,m^) V 

Using this result, the second amplitude can be simplified to read 


2ml 
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(75) 
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V(m| -m^)V 


(76) 


while the decay amplitudes for the processes Z\\ —)■ and Z^ —)■ U()'hF|| , in which the 

final-state meson is transversely polarized, are suppressed relative to fl76|l by factors of mvjmz 
and mv/mw, respectively. These power-suppressed amplitudes contribute to the decay rate 
at 0{ml/ml^) and are thus negligible. They will be neglected below. Up to very small 
corrections, the two amplitudes in 0741) thus have an identical structure. 

Similar to fl55|) . we can write the most general form of the decay amplitudes in the form 


iA{Z M+1U-) = ± 


g^fM 


4\/2 cos 9 


w 
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nir 
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k ■ q 
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z ly pM 
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(77) 


The last structure inside the parenthesis contributes only if the W and Z bosons are longitu¬ 
dinally polarized. It has no analog in the case of radiative Z or IF decays. 

Summing (averaging) over the polarizations of the final-state (initial-state) particles, and 
setting the meson mass to zero, we obtain from 0771) the decay rate 


r(Z ^ M+IU") = 


a{mz) mzf^ 
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Notice the close similarity of this result with fISB]) . The differences are that we must now 
evaluate the coupling a at the electroweak scale, where for simplicity we do not differentiate 
between mz and mw- Also, a phase-space suppression factor arises, and the third form factor 
yields a contribution that is absent in (l56|) . In deriving the above relation we have used 
that = ^! sin^ Oy/ and sin^ 9w = ^ — "rnwlrn\. 

Evaluating the first two Feynman graphs in Figure [HI we obtain 


= / dx(t)M{x,ii) 

Jo 

F^^ = 2 + [ dx(f)M{x,fi) 


+ 


X + {1 — x)r {1 — x) + xr 

Zd Zu 


-^ 3 ^= / dx fi) {I — 2x) 

Jo 


X + {1 — x)r {1 — x) + xr 
Zd Zr, 


(79) 


+ 


a; -|- (1 — x) r (1 — x) -|- xr 


Here r = rri^lm% = cos^ 9w, and Zg = Vg + ag are the left-handed couplings of quarks to the 
Z boson. Explicitly, we have Zu = J — |sin^ 6 'vi/ and Zd = —\ + |sin^ 6 *vi/. From the results 
shown in and fl76l) . we see that the local contribution from the third diagram in Figure |9] 
adds the term 2 to the form factors F^. Note that the relevant combination 


FIfM 

2 “ 


1 — r 




1 


(80) 


entering the total decay rate is independent of the form of the LCDA. It follows that the 
third (longitudinal) term in the expression fl78|l for the decay rate yields l/(2r). The integrals 
over the LCDA in our expressions for the form factors Ff\ can be evaluated analytically to 
any given order in the Gegenbauer expansion, but in practice it is easier to evaluate them 
numerically. Because the functions in the denominators are slowly varying with x, we find 
that higher-order Gegenbauer moments in the expansion flT^ of the LGDA fi) give very 

small contributions. This is the essence of the approach by Manohar 2^, which we illustrate 
in Appendix [Dl We obtain 


+ — = 1.911 + 0.003 af(/i) - O.Ollaf (/x) + 0(10-^), (81) 

I I I I 2r 

where the leading term 1.911 can be written in closed form as 

1 1 — “t“ 2t In T 

(l - 2r-MOr^)/^(r) - 12 r/(r)4-F — , with /(r) =- - ---. (82) 

\ / 2r (1 — r)'^ 

The calculation of 0{as) corrections to the results in fl79D is an interesting project, which we 
leave for future work. We note, however, that in this case the value of the factorization scale 
should be taken lower than the mass of the decaying Z boson. Two natural choices would 
be the typical momentum transfer, /i^ ~ 2fc ■ g = (m| — m^) ~ (43GeV)^, and (twice) the 
energy of the hnal-state meson in the Z-boson rest frame, fi ~ (m| — rn^)/mz ~ 20 GeV. In 
our analysis we shall use latter value, but since in the tree-level calculation presented here the 
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Decay mode 

Branching ratio 

ZO ^ TT^W^ 
ZO ^ pippT 

ZO ^ K^W^ 
ZO ^ K*^W^ 
Z^ D,W^ 
ZO ^ D^W^ 

(1.51 ±0.005/) ■ 10-1° 

(4.00 ±0.15/) ■ 10-1° 

(1.16 ±0.01/) ■ 10-11 
(1.96 ±0.12/) ■ 10-11 
(6.04 ± 0.20ckm ± 0.22/) ■ 10-i° 
(1.99 ± O.Mckm ± 0.10/) ■ 10-11 


Table 8: Predicted branching fractions for various Z —>■ MW decays, including error 
estimates due to the uncertainties in the CKM matrix elements and the meson decay 
constants. Uncertainties in the shapes of the LCDAs have a negligible impact at tree 
level. Not shown are perturbative uncertainties due to the neglect of 0{as) corrections. 


scale dependence enters only through the tiny corrections involving af^ 2 (h) 118ip this choice 

has no noticeable effect. 

Our predictions for the various Z —)■ MW branching ratios are collected in Table [HI They 
are smaller than the corresponding branching fractions for W —)■ My decays by more than a 
factor 20 for light mesons and by more than a factor 60 for heavy mesons. Note, however, the 
curious fact that the Z^ —)■ branching ratio is about 15 times larger than the Z^ —)■ vr^y 

branching ratio. The latter quantity is tiny, because it is proportional to ~ 7 ■ 10“^. 
We do not quote uncertainties related to the shapes of the meson LCDAs, which have a 
negligible impact on our result, see flHTp . We stress that to the quoted errors one should add 
an uncertainty accounting for our neglect of QCD radiative corrections. We expect that they 
can change the branching ratios by 10-20%. 


It is interesting to compare our results with those obtained by Manohar in j25|]. To this 
end, we used flHop to eliminate the parameter v and expand the function given in fl82p in powers 
of the weak mixing angle. This yields (with = sin^ 9w) 


r{z M+W-) = 
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ASmz 
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+ 4 +0.003 af + 


(83) 


where we do not show terms of 0{s^) and contributions from higher Gegenbauer moments. 
The leading contributions in this expression are in agreement with those found in [^, where 
an anlogous formula with the parenthesis replaced by ~ f + I + I ^tv + • • • is given. 


5 Experimental considerations 

Having obtained detailed and accurate theoretical predictions for a large class of very rare, 
exclusive radiative decays of Z and W bosons, we now address the question of how to search 
for such decays experimentally. While we do not perform an exhaustive feasibility study here, 
we discuss several ideas related to possible experimental analysis strategies. The goal is to get 
a feeling for how difficult it will be to observe some of the decay modes discussed in this work at 
present and future particle accelerators, and what accuracy one will be able to reach. We hrst 
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consider the LHC, where about 10^^ Z bosons and 5 ■ 10^^ W bosons will have been produced 
in both ATLAS and CMS by the time the high-luminosity run with an anticipated integrated 
luminosity of 3000 fb“^ has been completed. We also consider a future lepton collider such as 


TLEP [l6|, where one can hope to produce about 10^^ Z bosons per year with a dedicated 
run on the Z pole, while 10^ W-boson pairs per year would be produced in a run at the WW 
threshold. A run just above the tt threshold would also produce very large WW samples. 

At the LHC one needs to worry about triggering and reconstruction. The trigger for the 
decays Z —)■ My can be based on photons and muons. The energy of the hnal-state photon is 
comparable to the energy of the photons produced in the Higgs-boson decay h —)■ 77 , where 
it has been demonstrated that such events can be triggered on. Muons, on the other hand, 
are produced only in some cases, in particular in the decay modes containing a vector meson. 
Reconstructing these event appears to be challenging but not impossible. Probably the most 
promising modes are Z <//'07 ^md Z —)■ T(nS ')7 followed by a fully leptonic decay of 
the heavy quarkonium state. The corresponding rates, however, are very small. If only the 
muon channel can be used, then the combined Z —)■ H 7 —)■ /i+/i “7 branching fractions are 
of order 5 • 10“® for J/'^ and 1.5 • 10“® for T(IS'). Thus, we can expect to trigger on several 
hundred Z —)■ J/'^y —)■ /r^/r “7 events and up to one hundred Z —)■ T( 1 S ')7 —)■ /r^/i “7 
events at the LHC. While there is no signihcant physics background we can think of, the 
combinatorial background may be substantial. Given these challenges, it is encouraging that 
ATLAS has recently reported hrst upper bounds (at 95% CL) on the branching fractions 
Br(Z ^ J/i/j 7 ) < 2.6 ■ 10 - 6 , Br(Z ^ T(l^) 7 ) < 3.4 • lO"®, Br(Z ^ T(25) 7 ) < 6.5 • lO"®, 
and Br(Z —)■ T(3S') 7 ) < 5.4-lO”® j^. Further dedicated experimental studies of these decays 
would be worthwhile. The yield of Z —)■ Ty events can be enhanced by about a factor 2 if one 
combines the T(nS') channels with n = 1 , 2 , 3, which have similar leptonic branching fractions 
(ranging between 1.9% and 2.5%). The case Z T( 4 S ')7 may be particularly interesting, 
since the T(4S') resonance can decay to a pair of B mesons, which gives rise to displaced 
vertices. It might be possible to achieve a larger effective rate for this decay mode by using 
highly efficient 6 -tagging methods. Observing Z —)■ My decays into other hnal states seems 
to be difficult at the LHC. Ideas for reconstructing highly energetic cj), p and u mesons have 
been presented in [l^ in the context of a study of exclusive h —)■ Hy decays. For light mesons 
decaying into two photons, such as and 7 , it might be possible to tell that there are more 
than one photon in the hnal state provided one of the photons is converted into an e^e” pair. 

The situation with W —?• My decays at the LHC seems less promising. In this case the 
hnal state contains a charged hadron, and it is not clear to us how one could reconstruct it in 
the high-multiplicity environment of a hadron collider. Perhaps the most promising case is the 
decay W —)■ hlgy, of which over 10,000 events should be produced. The problem is how to tag 
the Dg mesons. A very interesting, dedicated study of several exclusive radiative W —?• My 
decays has been performed in [l^, to which we refer the reader for more details. 

It looks much more promising to search for rare exclusive Z My decays at a future 
lepton collider such as TLEP, in particular if one envisions a dedicated high-luminosity run on 
the Z pole. The advantage of such a “Z factory” is that it would produce yields of up to 10^^ 
Z bosons per year in a very clean environment. The detectors of such a future experiment 
are expected to have excellent particle-identihcation systems. This would make it possible to 
perform precision studies of many of the decay modes we have discussed. In particular, all 
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Z —)■ M'j decays except for Z 7r°7 should be accessible, and in several cases it should be 
possible to measure the branching ratios at the percent level. There is even hope for observing 
some of the weak radiative decays Z —)■ MW, whose branching fractions can be of order 
several times 10“^°. With dedicated runs above the WW or ti thresholds one would produce 
samples of W bosons large enough to search for the Cabibbo-allowed W —?• My decay modes. 
It would be most rewarding to perform detailed feasibility studies for measurements of rare 
exclusive Z and W decays at such a facility. 


6 Summary and conclusions 


Based on the formalism of QCD factorization, we have performed a comprehensive and sys¬ 
tematic analysis of the very rare, exclusive radiative decays of Z and W bosons into final 
states containing a single pseudoscalar or vector meson. The basis of our study is a factoriza¬ 
tion theorem derived in soft-collinear effective theory, which expresses the decay amplitudes 
as convolutions of calculable hard-scattering kernels with light-cone distribution amplitudes 
(LCDAs), in a systematic expansion in powers of (Aqcd/’^z,w)^ and {ruMl'i^z,wY■ For the 
hrst time, we have included the complete set of one-loop QCD radiative corrections. Large 
logarithms involving the ratio of the electroweak scale to the typical hadronic scale po, at 
which model predictions for the LCDAs are obtained, are resummed to all orders of perturba¬ 
tion theory. We have also estimated the leading power-suppressed effects, finding that their 
impact on the branching ratios is typically of order 10“^ compared with the leading terms. 
Larger power corrections up of to 1% can only arise for mesons containing heavy b quarks. 
The exclusive decays Z —)■ My and W —)■ My therefore offer an ideal laboratory for testing 
the QCD factorization approach in a controlled and theoretically clean way. Our main phe¬ 
nomenological results for the relevant branching ratios are collected in Table O We have not 
considered in this work the interesting decays Z —)■ y^'^y. Their analysis is complicated by 
the fact that the rj and rj' mesons contain a flavor-singlet Fock component, and to predict the 
decay rates at 0(as) one needs to take into account the two-gluon LCDA of these mesons. 


This will be discussed in a separate publication [33 


Our results form the basis of a rich, novel program of electroweak precision studies of 
the Standard Model, which offers powerful tests of the QCD factorization approach in a 
situation where it should deliver precise predictions. With the statistics obtainable in the 
high-luminosity run at the LHC and, in a much cleaner environment, at future lepton colliders, 
it will be possible to observe several of the very rare decays discussed here and measure their 
branching fractions with some accuracy. Precise rate measurements, which would be possible 
at a dedicated Z-boson factory, would offer the unique possibility to extract highly non-trivial 
information about the LCDAs of various mesons in a completely model-independent way. 
More specifically, for each meson M one will be able to extract the sums over the even and 
odd Gegenbauer moments, ^Fe electroweak scale fi ~ mz, 

up to small and calculable radiative corrections. This is a consequence of the structure of 
the basic convolution integrals in 0391) . We cannot imagine a theoretically cleaner way to get 
access to this kind of information. We have also performed an exploratory study of the weak 
radiative decays Z —)■ MW, which allow for tests of the QCD factorization approach at lower 
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Decay mode 

Branching ratio 

Decay mode 

Branching ratio 

Z^ —)■ 7r°7 

ZO ^ pOy 

ZO ^^7 

ZO ^ 07 

ZO ^ J/07 
ZO ^ T(1R)7 
ZO ^ T(4R)7 

(9.80 ± 1.03) • 10-^2 
(4.19 ±0.47) ■ 10-9 
(2.82 ±0.41) ■ 10-® 
(1.04 ±0.12) ■ 10-® 
(8.02 ±0.45) ■ 10-^ 
(5.39 ±0.16) ■ 10-^ 
(1.22 ±0.13) • 10-^ 

—)■ 7r=*=7 
IU± ^ p±7 
IU± ^ A±7 

^ A*±7 

1U± ^ D,7 
1U± ^ D±7 
1U± ^ R±7 

(4.00 ±0.83) ■ 10-9 
(8.74 ± 1.91) ■ 10-9 
(3.25 ±0.69) ■ 10-19 
(4.78 ± 1.15) ■ 10-19 

(3.66;if).iQ-8 
(1.38 ^9Ji). 10-9 
(1.55 ;9;79). 10-12 


Table 9: Summary table of our predictions for the branching fractions of exclusive 
radiative decays of Z and W bosons. Different sources of theoretical errors have been 
added in quadrature. 


scales /i ~ 10 GeV, which are only a few times higher than those relevant to exclusive hadronic 
i?-meson decays. Our predictions for the corresponding branching fractions obtained at tree 
level have been given in Table El 

Several generalizations and extensions of our work are possible and worth exploring. Our 
formalism can be applied in a straightforward way to obtain high-precision predictions for 
exclusive radiative (and weak radiative) decays of the Higgs boson, extending previous tree- 
level analyses presented in [l7-22|. One goal of such studies is to search for enhanced Yukawa 
couplings and flavor-changing interactions of the Higgs boson. In this context, new-physics 
studies analogous to those presented in Section 13.41 are particularly interesting. Without 
further conceptual developments, our formalism can also be extended to calculate the rates 
for purely hadronic decays, such as Z,W,h —)■ Mi M 2 or even decays with more than two 
particles in the hnal state. These extensions are left for future work. 

The physics case for studying some of the very rare, exclusive decays of heavy electroweak 
bosons is compelling to us. There is some beautiful physics to be explored here, both from 
the theoretical and the experimental points of view. We hope that our detailed exploratory 
survey will raise sufficient interest that some dedicated feasibility studies for discovering such 
decays at the high-luminosity LHC and future lepton colliders will be performed. 
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A Light-cone projectors for vector mesons 

The LCDAs of vector mesons at leading and subleading twist have been studied in great 
detail in The corresponding momentum-space projectors were derived in ji^. In 

analogy with ([S]), the light-cone projector at leading and subleading power for a longitudinally 
polarized vector meson reads 




Mv^i (fc, x,/i) = — —^(j)v{x,^i)- ^ 
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In the approximation where three-particle LCDAs are neglected, the QCD equations of motion 


imply the relations |4ll. |42 

/i|| (x, p.) 1^^^ (2x 1) 4’^(x, p), h|| (x, p) 1^^^ 24’^,(x, p), 


iy ^v(.y,y) - =i(i, 
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WWA 
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1-2/ 
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(A.3) 


In this approximation, the twist-3 two-particle amplitudes can be expressed in terms of the 
twist-2 LCDA (j)y. 

The light-cone projector for a transversely polarized vector meson is yet more complicated. 
Up to twist-3 order, one obtains 


MvAk,x,y) = !^^v* 0v(x,p) - ^ 
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Note that, compared with j^, the terms multiplying the Levi-Civita tensor in the second line 
have the opposite sign, because we are using a different sign convention for this object. In the 
approximation where three-particle LCDAs are neglected, the equations of motion yield the 
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relations 40-42 
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In this approximation, the twist-3 two-particle amplitndes can be expressed in terms of the 
twist-2 LCDA (j)y. 

As an application of these resnlts, we present the general expressions for the form factors 
entering the Z —)■ V±'j and W~^ —)■ V^'-^ decay amplitndes in fl45|l and ([59]). In the hrst 
case, we obtain (with x = 1 — x, and ignoring qnark-mass effects for simplicity) 
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We omit the scale dependence of the various quantities for simplicity. In the case of (15^ we 
must omit the prefactors Qy /Q and instead assign charge factors Qu and Qd to the two terms 
under each integral over x. These results can be simplihed signihcantly by using the relations 
flA.5p . The hnal expressions have been given in flT7|l and fl6T|l . 


B Determinations of meson decay constants 

We follow and determine the relevant meson decay constants from experimental data. The 
decay constants of charged pseudoscalar mesons can be determined from their semileptonic 
decays P~ —?• This analysis is performed by the Particle Data Group ji^, and it leads 

to the values for /p, fjy and f^s shown in Tables [D and [2J 


39 











































The decay constants of light charged mesons can also be obtained from the one-prong 
hadronic decays of the r lepton. The corresponding decay rates are given by 


r(r- ^ M-Ur) = 


EW 
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1 - 
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Vij are the relevant 


where 6p = 0 for pseudoscalar mesons and by = 2 for vector mesons 
CKM matrix elements. The factor S'ew = 1-0154 includes the leading-logarithmic |94J. 195 
and non-logarithmic electroweak corrections j^. From the measured branching fractions 
Br(7r-) = (10.83±0.06)%, Br(iF-) = (0.70±0.01)%, Br(p-) = (25.22±0.33)% and = 


(1.20 ± 0.07)%, along with the r-lepton lifetime r^- = (290.3 ± 0.5) -10 ^®s |45|, we extract 
U = (130.3 ± 0.4) MeV, fx = (154.3 ± 1.1) MeV, fp = (207.8 ± 1.4) MeV, and /k* = (203.2 ± 
5.9) MeV. The values of and fx are in excellent agreement with the (more precise) values 
extracted from semileptonic decays, supporting the reliability of this method. 

The decay constants of neutral vector mesons can be extracted from their electromagnetic 
decay width using 


r(V°-)• e+e ) = r(V°-)■/i+/i ) = a^{mv) cy , (B.2) 

oUly 

where the coefficients cy = related to the electric charges of the quarks 

that make up the vector meson [^. Explicitly, one has Cp = 1/2, Cuj = 1/18, = 1/9, 

cj/ip = 4/9, and cr = 1/9. The dominant QED corrections are accounted for by using the 
electromagnetic coupling evaluated at /r = my, which we compute using a(mz)~^ = 127.94 
and the approach described in [^. Averaging over the and yiJ,~ modes [4^, we obtain 
the measured branching fractions Br(p°) = (4.715±0.049)-10“^, Br(a;) = (7.284±0.140)-10“^, 
Br(0) = (2.952 ± 0.030) ■ lO"'^, and when combined with the total widths r(p°) = (147.8 ± 
0.9) MeV, r(a;) = (8.49T0.08) MeV, r(0) = (4.266T0.031) MeV this yields the decay constants 
fp = (216.3 ± 1.3) MeV, f^i = (194.2 ± 2.1) MeV and = (223.0 ± 1.4) MeV. For the heavy 
quarkonium states it is advantageous to use the measured electromagnetic width directly. 
They are Fee(J/V') = (5.55 ± 0.14) keV, Fee(T(lA)) = (1.340 ± 0.018) keV, Fee(T(2A)) = 
j^l2 ± 0.011) keV, Fee(T(35)) = (0.443 ± 0.008) keV and Fee(T(45)) = (0.272 ± 0.029) keV 
jidj ]. This gives the decay constants /j/^ = (403.3 ± 5.1) MeV, fr{is) = (684.4 ± 4.6) MeV, 
/t( 2 S) = (475.8 ± 4.3) MeV, /x( 3 S) = (411.3 ± 3.7) MeV and /t( 4 S) = (325.7 ± 17.4) MeV. 
For our analysis the value {Yln=i ~ ^ 4) MeV is also needed. As discussed 

in [^3 the combined effects of p-oj and u-(j) mixing lower by about 9.5 MeV and raise 
by about 7.6 MeV, while they have a negligible impact on fp. Note also that there are 
some uncertainties related to the sizable width of the p resonance, which we have ignored 
here. They might explain why the value of fpo extracted here is larger than the value of fp- 
extracted above from r decay. Combining the various extractions and using conservative error 
estimates, we obtain the values shown in Table [U 

^We are grateful to R. Zwicky for pointing out a numerical mistake in the published version of this paper. 
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C Gegenbauer expansion of the convolution integrals 


In order to derive closed analytic expressions for the basic convolution integrals and 
defined in (15^ . we employ the definition of the Gegenbauer polynomials in terms of the 
generating function 
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(C.l) 


From (IT^ and the second relation in (15^ . it follows that the coefficients defined 
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(C.3) 


The coefficients in fl4n|l are the coefficients of t"^ in the series expansion of h{t) 

around f = 0. It is then a straightforward exercise to derive expression flTTl) . 


D Connection with the approach by Manohar 

The approach put forward by Manohar in was originally developed as a tool to study the 
exclusive decay Z —)■ ITtf by means of a local operator-product expansion. It corresponds to 
the following series expansion of the propagator in the hrst diagram in Figure O 
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(D.2) 

oz + 

Using this expansion, the convolution integrals over the pion LCDA in flTOjl can be expressed 
in terms of local operator matrix elements, and the leading terms are determined by the 
normalization of the LCDA. For the phenomenological estimates in 25| and 15| only the 
leading term was kept. 

When the same approach is used for radiative decays, the IF-boson mass in the above 
relations must be set to zero, in which case one obtains 
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a:m| 
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m| 


1 + 5^(1-2x)” 
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(D.3) 
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Keeping the leading term is now unjustified, but if it is done this corresponds to replacing 
l/x —)■ 2 under the convolution integrals in ([39]). In the asymptotic limit, where all Gegenbauer 
moments are set to zero, these integrals are therefore too small by a factor 2/3. This explains 
the discrepancies in the predictions for the Z —)■ Mtt and hh —>■ Mtt branching fractions 
mentioned in Section IXTl 
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